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Summary 


Simulation of high speed propulsion systems may be divided into two categories, nonlinear and linear. The 
nonlinear simulations are usually based on multidimensional computational fluid dynamics (CFD) method- 
ologies and tend to provide high resolution results that show the fine detail of the flow. Consequently, these 
simulations are large, numerically intensive, and run much slower than real-time. The linear simulations 
are usually based on large lumping techniques that are linearized about a steady-state operating condition. 
These simplistic models often run at or near real-time but do not always capture the detailed dynamics of the 
plant. Under a grant sponsored by the NASA Lewis Research Center, Cleveland, Ohio, a new method has 
been developed that can be used to generate improved linear models for control design from multidimensional 
steady-state CFD results. This CFD-based linear modeling technique provides a small perturbation model 
that can be used for control applications and real-time simulations. It is important to note the utility of the 
modeling procedure; all that is needed to obtain a linear model of the propulsion system is the geometry and 
steady-state operating conditions from a multidimensional CFD simulation or experiment. This research 
represents a beginning step in establishing a bridge between the controls discipline and the CFD discipline 
so that the control engineer is able to effectively use multidimensional CFD results in control system design 
and analysis. 

1 Introduction 

The development of inlet models for high speed propulsion systems is important because of the current 
interest in high speed air-breathing propulsion systems. Modeling of these systems is difficult, because the 
complex physical processes are represented by nonlinear partial differential equations (PDE). An accurate 
plant model is required to develop a control system for the plant; the more accurate the model, the better 
the control design. Typically these models are either based on traditional propulsion control models or CFD 
models. 


Traditional propulsion control models typically utilize a large lumping technique for the spatial derivatives 
so that the propulsion system is represented by a set of nonlinear ordinary differential equations (ODE). 
These equations are often linearized about a steady-state point so that the control model is linear. Methods 
based on this linear ODE approach have been developed for propulsion systems, some of which are: the Cole- 
Willoh model [1], the Martin model [2], the Barry models [3], circuit models [4], and the Laplace transform 
of the Green’s function method [5] [6]. Unfortunately these models are often difficult to implement, do not 
always capture the nonlinear dynamics of the system, and are not typically used for multidimensional flows. 

Accurate nonlinear models of complex flows are usually obtained from CFD codes [7] [8]. These models 
can to some degree predict the behavior of large perturbations in the flow field, including unstart, buzz, 
turbulence, boundary layer growth, et cetera. Typically these CFD models are based on a large number of 
nodes which can then be used in a finite difference method to produce a large system of nonlinear equations. 
However due to their nonlinearity and large size, these models require large amounts of computational time 
and therefore are not suitable for controls analysis and design. 

An effective propulsion system model for control system design must adequately capture the dynamics 
of the system but also be of small order. CFD models fulfill the first requirement, and traditional controls 
models fulfill the second requirement. Therefore, a method that is based on both ideas might provide a 
reasonable model for controls applications. This concept has already been illustrated for one dimensional 
CFD models [9], In that paper, the development of a CFD-based linear modeling method combined with 
model reduction is used to model the inviscid flow of an axisymmetric one dimensional mixed compression 
inlet. 

In this paper, the CFD-based linear modeling method is applied to the inviscid flow of an axisymmetric 
multidimensional mixed compression inlet. It should be noted that the whole inlet needs to be modeled in 
order to accurately capture transient behavior. The CFD code PARC is used to obtain the steady-state and 
transient data. The inflow boundary condition is assumed to be supersonic, the outflow boundary condition 
is, assumed to be subsonic, and the Mach number at the exit is used as the boundary condition input. 
The next section describes the multidimensional CFD model development which is the basis for the linear 
model of the inlet; this includes the governing equations, the development of the split flux model, and the 
development of the boundary conditions. Then in section three, the linear model is derived by implementing 
linearized methods of the previous section; an input and various outputs are also developed in this section. 
A summary of the model reduction method and calculation of the associated error bounds follow in sections 
four and five. In section six, an application of the method is illustrated on a mixed compression inlet, and 
a conclusion follows in section seven. 


2 CFD Model Development 


2.1 Governing Equations 


The dynamics of an internal flow propulsion system are often represented by the nonviscous multidimen- 
sional Euler equations. The conservative form of these equations is defined by Hirsch [10] as: 


Conservation of Mass: 

d(p) , d(pu) , d(H d{pw) 

dt dx dy dz 


( 2 . 1 ) 


Conservation of Momentum: 

d{pu ) d ( pu 2 + p) d ( puv ) d ( puw ) 

dt dx dy dz 


2 


0 


( 2 . 2 ) 


d ( pv ) d ( puv ) d ( pv 2 + p) d (pvw) 

dt dx dy dz 

d ( pw ) d ( puw ) d ( pvw ) d (pw 2 + p) 

dt dx dy dz 


0 


Conservation of Energy: 

djpE) d[u(pE + p)\ i d[v(pE-\-p)\ d[w(pE + p)\ 


dt 


+ 


dx 


+ 


dy 


+ 


dz 


= o 


( 2 . 3 ) 


The number of equations in the nonlinear system of partial differential equations is SiViA^A^ where 
and Ns are the number of grid points in the x, y, and ^-directions. The two dimensional form of the Euler 
equations is obtained by removing the third momentum equation and the spatial derivatives with respect to 
z in the remaining equations. Throughout this paper the three dimensional form of the equations will be 
implemented. 

If the partial derivative terms are expanded, the conservative form of the equations can be rewritten as: 

dli df (Tl) d~g (~u) dli (It) 


+ 


+ 


+ 


= 0 


dt dx dy ' dz 

where the vector components, it, f (it) , if (it) and h (It) , are defined below. 
State vector: 


u = 


( 2 . 4 ) 


p 


P 

pu 


mi 

pv 

= 

rn 2 

pw 


ms 

pE 


£ 


The total energy per unit volume, e, is defined as, £ = pE — p ( e + 
energy per unit mass. 


u z -f v z + w z 


( 2 . 5 ) 


, where e is the internal 


Flux vectors: 


f C*?) = 


pu 

pv? +p 
puv 
puw 

L u(e+p) J 


( 2 . 6 ) 


~9 (~u) = 


pv 

puv 

pv 2 +p 
pvw 

v(e-\-p) 


( 2 . 7 ) 
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pw 
puw 

h ( ~u ) — | pvw 
pw 2 + p 

L ™(e+?) J 

For a perfect gas the static pressure can be defined as: 


( 2 . 8 ) 


V = (7~l)(e-^ 2 ) 

= (t-1) 


= (7 - 1) f e~p- 


u 2 + v 2 4- u> 2 


(2.9) 


where !? = ui + u j + wk, and the flux vectors may be rewritten in terms of the state variables. Since the 
flux vectors are homogenous functions of degree one in ~u } they can be written as: 

/ (u) = A~u 


~g {li) = BJl 


(2.10) 


h ( 111 ) — Cli 

Now the derivative of the first flux vector with respect to x can be represented in quasi-linear form by the 
following, 


where A is the Jacobian of the flux vector. 


d f (~u) _ , dli 
dx dx 


df (IT) 

dli 


( 2 . 11 ) 


( 2 . 12 ) 


The Jacobians of the other flux vectors may be computed in the same manner. The substitution of the 
Jacobians into equation (2.4) results in the following partial differential equation, 


dli .dli n d~u ~,d~u 

"or + A—~ + B~~ + C— = 0 
dt dx dy dz 


(2.13) 


where the Jacobians are defined as: 


2 , 2 1-^2 


0 


0 


0 


H — v 


(3 — 7) u 


—uv 


—uw 


(7 — 1) v — (7 — 1) w 7 — 1 
u 0 0 
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0 
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2? - (7 - 1 ) V 2 
p 


ye 7 — 1 


(l ? 2 + 2u 2 ) — (7 — 1 )uv — (7 — 1) uw 7 u 


(2.14) 
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B = 


~uv 
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-V 2 + ~ o ~ ~^ 2 — ( 7 — l)u 


c = 


2 

-vw 


L f-(7-DV 2 


~uw 


_^ 2 + l_J:V 2 


(3 - 7) v 
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0 0 

— (7 — 1 ) w 7 — 1 
v 0 


(7 — 1 )uv — — (l? 2 -f 2 u 2 ) — (7 — 1 ) vw jv 

pi 


0 


0 

0 
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(7 — 1 ) u — (7 — 1 ) v 


(3 — 7 ) w 


0 

0 

0 

(7- 1) 


7£ 

L p 


- (7- 1 ) 


y£ y — X 

— ( 7 — l)^ — (7 — 1 )vw — ( ~v* 2 + 2w 2 ) 7 w 

P ^ 


The characteristics, or local eigenvalues, of the first Jacobian are equal to: 


the characteristics of the second Jacobian are equal to: 


and the characteristics of the third Jacobian are equal to 


^1,2,3 

- 

u 

a 4 

= 

u + c 

A 5 

= 

u — c 

e equal 

to: 


Al,2,3 

= 

V 

a 4 

— 

V + c 

As 

= 

v — c 

are equal 

to: 

Al,2,3 

= 

VO 

a 4 

= 

w 4 - c 

As 


w — c 


(2.15) 


(2.16) 


(2.17) 


(2.18) 


(2.19) 


Note that c is the speed of sound. A system of ordinary differential equations that approximate equation 
(2.13) can be obtained by replacing the spatial derivative terms with finite difference expressions; then 
the system of equations may be integrated numerically to obtain the flow field solution. In order for the 
overall system to be numerically stable, the direction of the characteristics must be taken into account 
when the spatial derivatives are replaced. For example in the axial direction when the flow is supersonic, 

the characteristics are all positive, and one finite difference expression can be used for the — — term. If 

the flow is subsonic, the signs of the characteristics are mixed, and a single finite difference expression for 
d V/ 

~q~ will create an unstable set of ordinary differential equations. If the Jacobians of equation (2.13) are 
split according to the signs of the characteristics, then different finite difference expressions for the spatial 
derivatives can be used for each of the positive and negative terms. The next section illustrates how to split 
the system into its positive and negative parts. 
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2.2 Split Flux Model 


The split flux method detailed in references [10] and [11] is summarized in this section. The split flux 
method separates a flux vector into subvectors which correspond to the positive and negative characteristics 
of the Jacobian. The split flux model can be written as the following equation, 


d'u df + (~u*) df (Ti) dg+(~u ) dg (Tl) dh + (Tl) dh (Tl) 
dt dx dx dy dy dz dz 

and the positive and negative subvectors can be calculated from the following: 

f^(Tt) = A±T? 


( 2 . 20 ) 


(Tl) = B ± Tl 


( 2 . 21 ) 


h^Cvt) = C^Tl 

Substitution of equation (2.21) into the split flux model equation (2.20) produces the following result. 


* +A+^+A-^+B+^+B-^ + &^+Cr^-= 0 

at ox ox oy dy oz oz 


The positive and negative Jacobians satisfy, 


and are calculated from: 


A = A+ + A- 
B =. B+ + B- 

c = c + + c- 


i± - 


±rr - 1 


KiAfKi 


B ± - K 2 A 0 K 0 1 


( 2 . 22 ) 


(2.23) 


(2.24) 


CP = K 3 A3 

The right eigenvectors of A, B , and C are defined as: 


Kx = 


1 

0 

0 




2c 

2 c 

u 

0 

0 

£ {u+c) 

l>- c > 

V 

0 


pv 

pv 

p 

Y 

Yc 

w 

p 

0 

pw 

pw 



2c 

Yc 

1 ? 2 


... P 1 

( ~v 2 c 2 \ 

1 1 ™ i \ 

p f ~v 2 c 2 


2c 


7-1 


+ “ 2J 


7 - 1 


(2.25) 
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k 2 = 


0 

0 

0 

-p 

—pw 


u 


w 


U > 2 


P 

0 

0 

pu 


P_ 

2c 

pu 

2 c 


2 c 


(v 4 - c) 


pw 

2 c 


_P 

2 c 

pu 

2 c 

pw 

~2c 


P_ 

2c 


1 T 2 




7“ 1 


+ cu ^ 


#3 


0 

-p 

0 

0 


u 


tc 


_P 

2 c 

pu 

2 c 

pv 

2 c 


" 2 " + 7 


_P 

2 c 

/m 

2 c 

pv 

2 c 


— cv 


2c 


{w + c) 


2 c 


(w — c) 


pv —pu 


p_ 

2c 


+ 


7-1 


+ cw ) — 


+ 


and the matrices Af , A^, and A 3 can be defined as: 


Af - 


^(w±|u|) 0 


§(«±M) 


Uv±\v\) 0 


s( v± M) 


A? = 


i(w±H) 


i 0±M) 


0 

0 

0 

0 

0 

0 

\ (v±H) 
0 
0 

0 

0 

J(w±H) 

0 

0 


k (u + c ± \u + c|) 


7-1 

0 

0 

0 

0 


cw 


k(u-c±\u-c\) 


i (v + c i \v -f cj 


0 


h(v-c±\v-c\) 


0 0 

0 0 

0 0 

\ {w -f c± \w + cj) 0 

0 ^ (w — c ± \w — c|) 


(2.26) 


(2.27) 


(2.28) 


(2.29) 


(2.30) 


There are a variety of splittings that can be used for A. As long as the characteristics of A + and A~ satisfy 
A = A 4 " + A - , the splitting is valid. Once the system, equation (2.22) , is split into its positive and negative 
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Jacobians, a different finite difference expression can be used to approximate the spatial derivatives for each 
Jacobian. The spatial derivatives associated with the positive Jacobians are discretized with a backward 
difference operator: 


dli' 


i,j,k 


dx 

9 U i t j t k 
dy 

& U i,j,k 


u i,j,k u 

•Ei,j,k 3'i—l,j,k 

M i,j,k U j,j— l,fc 

Vi,j,k Vi,j—l,k 


(2.31) 


dz 


U i,j,k u i,j,k— 1 

z i,j,k ~ z i,j,k~ 1 

and the spatial derivatives associated with the negative Jacobians are discretized with a forward difference 
operator: 

^ i+l,j,k ^ i,j,k 


dli 


ijjyk 


dx 

® u i,j,k 
dy 

d U i t j t k 
dz 


x i+l,j,k x i,j,k 
U i,j+l,k ^ i>j,k 

yi,j+l,k yi,j,k 

H i,j,k + 1 U i,j,k 


(2.32) 


z i,j,k+l z i,j,k 

The grid point is denoted by the subscript z, j, k. The approximations for the spatial derivatives are substi- 
tuted into equation (2.22) which results in the following equation at each grid point of the system. 


0 = te_r:T )+^ 


x i,j,k %i—l,j,k 


V- i+l,j,k U i t j,k 
x i+l,j,k x i,j,k 


+ B Zj,k 


U i,j,k U i,j—l,k 

Vi,j,k yi,j—l,k 

— y — ► 

U iJ>k ^ i,j,k— 

z i,j,k z i,j,k — 1 


+ B- 


i,j,k 




^ ^ i,j,k 

yi,j+l,k yi,j,k 

u i,j,k + 1 u i,j,k 

z i,j,k - fl — z i,j,k 


(2.33) 


This can be rewritten as: 
9 U i,j,k 


0 = 


^i,j,k 


dt 


x i,j,k ^i—l,j,k 


+ 


A + 

A i,j,k 


^ i—l,j,k 


A i,j,k 


K i.k 


Vi,j,k yi,j~l,k 


U i,j—l,k 




z i,j,k z i,j,k— 1 


^ i,j,k—l 


+ 




x i,j,k Xi— l t j t k x i,j,k yi,j,k yi,j — l,k 


B i,j,k 


+ 


c Lk 


(2.34) 


^ i,j,k 


Vi,j+l,k yi,j,k z i,j,k z i,j,k— 1 z i,j,k + 1 z i,j,k 


U i,j>k 


+ ■ 


Ai,j,k 


^ i+l,j,k + 


B i,j,k 


V* i,j+l,k T 


^i,j,k 


U i,j,k+l 


x i+l,j,k x i,j,k yi,j+l,k yi,j,k ’ z i,j,k + 1 — z i,j,k 

Equation (2.34) represents the dynamics of the internal grid points of the CFD model; there are still boundary 
conditions that must be satisfied at the following locations: 

i = j — k — l 
i = Ni 

j = n 2 

k = N 3 


8 



where N\, N 2 , and iV 3 are the total number of grid points in the x, y and z-directions. The boundary 
conditions are developed in the next section. 

2.3 Boundary Conditions 

Boundary conditions can be categorized as either physical or numerical. Numerical boundary conditions 
correspond to characteristics leaving the domain; therefore they are determined from the interior grid points. 
The physical boundary conditions correspond to characteristics entering the domain and cannot be deter- 
mined from the interior grid points; therefore, they must be specified. The numerical treatment for the 
boundary conditions follows in the next two sections. 


2.3.1 Wall Boundary Conditions 


The boundary conditions at the y and z-planes are implemented using the method of non-reflective 
boundary conditions [10]. When using this method, the physical boundary conditions are set equal to zero, 
and the numerical boundary conditions are determined from the interior grid points of the computational 
grid. Since the characteristics at the boundary are propagating in one direction, one finite difference equation 
can be used to replace the spatial derivative. The following is a general equation that can be used at these 
boundaries (j = 1, j =■ N 2 , k = 1, and k — N 3 ) for the boundary conditions. 


0 = 


dlt 


i,j,k 


dt 


a + 

A i,3,k 

x i,j,k x i—l,j,k 


^ i— l,j,k "f - 


A + 

A i,j,k 


i.J.k 


x i,j,k %i— l,j,k %i-\r l,j,k x i,j,k 


u i,j,k 


+ 


Aj,j,k 


(2.35) 


, i+ „* - + + 


At the boundaries, A 2hc and A 3 bc are determined from the numerical boundary conditions; where as, A + and 
A~ are determined from the positive and negative Jacobians from the split flux method. 


d u 


For j — 1, the spatial derivative is replaced with a forward finite difference, equation (2.32) , hence 


d~u 


A 2bc must have all negative characteristics. Therefore, — , A 2bc , B + , and B~ become the following: 


dli 

dy 


^ i,j+l,k U i,j,k 

yi,j+l,k yi,j,k 


^2 bc 


0 0 0 0 o' 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 v — c 


v > 0 


A 


2 be 


V 0 0 0 0 

0 V 0 0 0 

0 0 V 0 0 

0 0 0 0 0 

0 0 0 0 v — c 


v < 0 


(2.36) 
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B + = 0 

B- = K 2 A 2bc K j- 1 


Likewise, for k — 1, 


— A r+ 

» A 3 bc . C , 


and C become the following: 


dH 5 

~dl 




A 3 6c = 


0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 w — c 


w > 0 


A 3 


be 


w 0 0 0 0 

0 w 0 0 0 

0 0 w 0 0 

0 0 0 0 0 

0 0 0 0 w — c 


w < 0 


(2.37) 


C+ = 0 

C- = K 3 A 3 ,J{^ 


d~u 

For j — N 2 , the spatial derivative is replaced with a backward finite difference, equation (2.31) , hence 

Q 'll __ 

A 2 bc must have all positive characteristics. Therefore, , A 2 bc , B + , and B~ become the following: 


H ^ i,j,k V* i,j—l,k 

Ui,j,k yi,j — l,k 


A2 h 


v 0 0 
0 v 0 
0 0 v 
0 0 0 
0 0 0 


0 0 

0 0 

0 0 

v + c 0 

0 0 


, v > 0 


a 2 


be 


0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 v + c 0 

0 0 0 0 0 


v < 0 


(2.38) 


B + = K 2 A 2hc K^ 

B~ = 0 
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Q 11 

and for k = N 3 , — — , A 3 bc , C + , and C~ become the following: 
dz 


dli 

Ite 


W i,j,k U 1 


z i,j,k z i,j,k — 1 


A 3 ()C = 


W 


0 w 
0 0 


0 0 
0 


w 


0 0 0 w + c 0 

0 0 0 0 0 


w > 0 


A-3 bc = 


0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 w + c 0 

0 0 0 0 0 


w < 0 


C+ - K s A 3bc K^ 

c~ = 0 


(2.39) 


The coefficients for the wall boundary condition equations listed in section 3.2 will be computed from the 
B~, C + , and C~ terms defined in the above equations. 

2.3.2 Inflow and Outflow Boundary Conditions 


Compatibility relations with time-differenced physical boundary conditions is the method used for the 
treatment of the inflow and outflow boundary conditions. Using this method, the incoming characteristics are 
replaced with specific physical boundary conditions so that only information transmitted from the interior 
is maintained. The general development for the treatment of the boundary conditions as taken from Hirsch 
[10] and Chakravarthy [12] follows. 


If the governing PDE, equation (2.4) from section 2.1, 


dli 

dt 


+ 


df (~u) 

dx 


d~9 {n) + d h (u) = Q 


dy 


dz 


is rewritten in terms of the characteristic variables for the time derivative, the characteristic variables can 
be split into the incoming physical characteristics, ur , and the outgoing numerical characteristics, w N . 


d_ 

dt 



+ 


a 

A—~ + 

' (xr'f ' 

B^l + 


w 


dx 

. (^r 1 ) P . 

dy 

_ ( K ^) P _ 


, dli 
dz 


— 0 


(2.40) 


The partitioning of the K\ matrix shown here will be discussed in detail in section 3.3.1. Note that the 
position of the physical and numerical characteristics will be different for different types of boundary condi- 
tions. With this boundary condition method, information transmitted by the characteristics to the interior is 
discarded. Therefore, the incoming characteristics are set equal to zero by replacing them with the physical 
boundary condition, Bb c = 0, as follows. 


d_ 

dt 


" ' 

+ 

‘ ^r 1 )" ' 

a 

A— (- 

' ( K i l f ' 

B if + 

" ( k ^) n ' 

B bc 


0 

dx 

0 

dy 

0 


C 


dli 

Ite 


= 0 (2.41) 
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Converting back to conservative variables and replacing with yields, 


(^rrf 

dBbc dli 
dli dt 


+ 


-1 \N 


( K D 

o 


OX 


dt 

-1 \N 


d u dt 


( K n 

0 


B^ + 

dy 


( K D 

0 


N 


c& i = °( 2 - 42 ) 


Now the implementation of the time-differenced boundary condition is performed. This is done by taking 
the derivative of Bbc with respect to time and replacing the derivative by a finite difference equation. 


dBbc 

dt 


dBbc dUt 
dli dt 


B i>c (V n+ 1 )-B 6 c (V’*) = (V n+1 - V") 

Bfc c (V n+1 ) = B 6 c (V") + |^(V’ !+ 1 - V”) 


Now Bbc (H >n+1 ) is set equal to zero, because it has replaced the incoming characteristic. 

0 = B tc (V") + ^(l? n+1 - V") 


-B bc {li n ) — 

Bbc ( 11 ) — 


dBbc 
dli 

dBbc dm 


(-^n+l _ 


du dt 

Substitute the last line of (2.44) into equation (2.42). 


Let L\ be defined as: 


and 1/2 be defined as: 


Li = 




(Ki 1 ) 

dB bc 


N 


d~u 

-1\N 


(*T ) 
o 


Substitute (2.46) and (2.47) into equation (2.45) and simphfy. 




Li 1 


0 

—Bbc 

0 

-Bbc 


(2.43) 


(2.44) 



dli 

( K lf 

,ent 

( k D n 

n d~u 

(Ki'f 

„d~u 

0 

1 

CD 

* o 

1 


0 


0 

5— + 
dy 

0 

c i^ = 

—Bbc 


(2.45) 


(2.46) 


(2.47) 


(2.48) 
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Since the incoming characteristics have been zeroed out and replaced with physical boundary conditions, the 
spatial derivative can be replaced with one finite difference expression. At the cowl lip the spatial derivative 
with respect to x is replaced with a forward difference equation (2.32), and at the compressor face it is 
replaced with a backward difference equation (2.31). For the spatial derivatives in the other directions, they 
are replaced after and have been defined either from the boundary conditions or from the split flux 
method, as shown in the following equations. 


Cowl lip boundary condition: 


dVl ' i ' k +LpL i 


dt 


A 1 t j t k — ► A\ t j t le — ► 

u 2,j,k ~ — u 1 ,j,k 


B U 


+ 


B L k 


^2 ,j,k 


B 1 ,j,k 


%2,j,k %l,j,k 


yi,j,k - yi,j-i,k 


u 1.J— l,fc 


yi,j,k - yi,j-i,k yi,j+i,k - yi,j,k 

✓ 

ri+ 

°1 ,j,k 


V* 1 ,j,k T 


B 1 ,3,k 


yi,j+i,k - yi,j,k 




z l ,j,k z l,j,k—l 


u + 


CTj,» 


z l,j,k z l,j,k—l z l,j,k+l Z 1 ,j,k 


u 1,3+1, k 


U 1, j,k 


(2.49) 


+- 


°Tj> 

z l,j,k+l ~ z l,j,k 


U ) “ -^l 


-1 


Bbc 

0 


Compressor face boundary condition: 


d~u 


N\,j,k r _i 


dt 


+ L^L 2 


A 


N u j,k 


x Ni,j,k x Ni—l,j,k 


^ Nl ,j,k 


A N u j, 


i,k 


x N\,j,k x Ni — l,j,k 


Ni—l,j,k 


R+ 

_Ni ,j,k 


yNj.,j,k - yN u j-i,k 


u Ni,j-l,k + 


B + 

n N lt j,k 


B 


Ni,j,k 


VN lt j,k ~ VNuj-l.k VNuj+hk ~ yN u j,k 


u Ni ,j,k 


B 


+- 


Ni ,j,k 


VNij+l,k ~ VN lt j,k 


u N lt j+l,k 


n+ 

° N u j,k 


z Ni,j,k z Ni,j,k— 1 


u Ni,j,k—1 + 


n+ 

,j,k 


Z Ni ,j,k z Ni,j,k - 1 


(2.50) 


C 


Nt ,j,k 


z Ni,j,k+l Z N\ ,j,k 


u N lt j,k + 


c 


Nr,j,k 


Z N\ ,j,k + 1 z Ni ,j,k 


u Ni,j,k+1 I — 


-1 


0 

-Bbc 


The next section develops the linear model using the split flux technique from section 2.2 and the boundary 
condition methods from section 2.3. 


3 Linear Model Development 

3.1 Small Perturbation Model 

For a small perturbation model, the states, inputs, and outputs of a nonlinear system are assumed to 
be the combination of a steady-state value and a small time dependent perturbation as shown below in 
equations (3.1) through (3.3). 

States: 

A* ~ ss + (3-1) 
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Inputs: 


V = Tf„ + sV 


(3.2) 


Outputs: _ 

3 ? =y„ + «y t ( 3 . 3 ) 

The small perturbation model is only valid when operating within a small region around the steady-state 
value; once outside this region, the linear model is no longer an accurate representation of the nonlinear 
system. In other words, the farther the nonlinear model moves away from the operating point, the less 
accurate the linear dynamic response. 


Because of the large number of equations in the CFD-based model, it is convenient to describe the small 
perturbation model in state space format; this also facilitates the placement of the inputs and outputs for 
the system. 

^-6lt = A81 + B8 V 

at 

sY = celt + D£E? (3.4) 

The A, B, C, and D matrices are defined as follows: 


A = 

system matrix 

B - 

input matrix 

C = 

output matrix 

D - 

input/output matrix 


The data needed for the calculation of these matrices is obtained from a steady-state operating point of the 
propulsion system model. The contents of these matrices are developed in the following sections. 


3.2 System Matrix 


The system matrix for the small perturbation model is generated by applying equation (2.34) to each 
grid point of the propulsion system with the Jacobians evaluated at the steady-state operating point. For 
the interior grid points, equation (2.34) can be rewritten as follows: 


d8 U i ! j l k 

dt 


— ^ i— l,j,k + b i,j,k8 U i t j—l 7 k “h Si,j,k8 U 1 T U i t j t k 

+Gi t j 7 kS U -f" C i t j t k8 U i 7 j-\-\ t k 4” f i,j,k8 U 


where the coefficients on 6~u are defined as: 

a: 


&i,j,k — 


i,3,k 


A + 

^i,j,k 


+ 


B i,j,k 


•Ei+l,j,k %i,j,k 3'i,j,k •^i—l,j,k Ui,j+i,k yi,j,k 




+ 


^ i,3,k 


c Lk 


H,j,k — 


C i>j,k — 


Vi,j,k yi,j—l,k Zi,j,k-\-l Zi,j,k %i,j,k ^i,j,k — 1 


yi,j,k yi,j—l,k 


B i,j,k 


yi,j+l,k yi,j,k 


(3.6) 
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( 3 . 7 ) 


&i,j,k ~ 

U,j,k — 
6 i,j,k — 


A + 
A i,3,k 


x i,j,k x i 

— 1 ,j,k 


,k 

x i+l,j,k ~ 

x i,j,k 


k 

z i,j,k+l ~ 

Z i,j,k 




z i,j,k z i,j,k— 1 

This can be written in the following compact form: 

d 


—SX = A SX 
dt 


( 3 . 8 ) 


If the i , j, and k indices are incremented in the following manner (((fc = 1, • • • , A/3) j = 1 , • ■ • , N2 ) i = 1 , • • • , N\) 
SX is defined as: 

S u 1,1,1 
^1, 1,2 

S^ 1 , 1 ,N 3 
S^ 1,2,1 




^^l,2,iV3 


^ 1,^,1 


( 3 . 9 ) 


S u i,n 2 ,n 3 

6 ~u 2,1,1 

. S~u n u n 2 ,n 3 . 

The banded structure of the A matrix is shown below: 

a f c e 

g a f c e 

g a f c e 

g a f c 

a _ b g a f c 

b g a f c 

b g a f 

d b g a f 

d b g a f 

d b g a 


( 3 . 10 ) 


The size of the system matrix is 5N1N2N3 by 5 Ah A/2 A/3. There may be some confusion between the system 
eigenvalues, which are the eigenvalues of the A matrix, and the local eigenvalues, which are the characteristics 
of the Jacobians. The local eigenvalues may be positive or negative depending on the flow conditions, but 
the eigenvalues of a stable continuous time linear model must have negative real parts. Therefore, all the 
eigenvalues of A must have negative real parts. The small perturbation model will be stable when the spatial 
derivatives of equation (2.22) are properly discretized according to the local eigenvalues. 
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The procedure for the development of the system matrix must be modified when the boundary conditions 
are taken into consideration. The linearized equations needed for the wall boundary conditions and inflow 
boundary conditions will be developed below. The linearized equations for the outflow boundary conditions 
are listed, but their development is detailed in section 3.3. 


3.2.1 Wall Boundary Conditions for the System Matrix 


The wall boundary conditions are implemented as modifications to the system matrix following the 
method developed in section 2.3.1. With i — 2, ■ • • ,Ni - 1, the equations and coefficients are shown below. 


d8~u 


i, i,i 


ft 

dt,i,i = 
^,1,1 = 

1,1 = 

f 4 ,1,1 = 


-1- c^i ,i6~Ui t 2,i + i,i,2 (3-11) 

A+ 

1,1 

x i, 1,1 x i— 1,1,1 

^z,l,l 


Bj, 1,1 _j_ Q, 1 ,1 


x i+ 1,1,1 — ^z.1,1 x i, 1,1 _ x i— 1,1,1 Vi, 2,1 — Vi, 1,1 Zi, 1,2 ~ Zi t \ t i 

~ A i,l,l 

x i- (-1,1,1 x i,l,l 


—B. 


*,i,i 


Vi, 2,1 — Vi, 1,1 
Zi, 1,2 -^i,l,l 


9S U 
dt 

di ( l ; fc 
6z,l,fc 
&d,l,k 
®z, l,k 
^i,l, k 
fi,l,fc 


di,l,k8~U i-l,l,k + gi,l,k8 ~Ui,l,k-l + &i,l,k8~U i,l,k (3.12) 

+ei,i,fc<5 u i-fi^fc + Ci,i t ic61i i t 2 ,k + f*,i,fc^'t?»,i l fc+i, k = 2, • • • , IV3 — 1 

A + 

A i, l,k 


x i,l,k x i—l,l,k 

z i,l,k ^z,l,fc— 1 
^z,l,fc 


Ki,i 


, -^z.l.fc 

_l l_ 


cri. fc c+, 


x i~\~ 1,1,^ x i,l,k x i,l,k x i—l,l,k Vi,2,k Vi,l,k ^Zjl^+l ^z,l,fc Zi^^ ^z,l,A:— 1 
~^z,l,fc 

^-z+l,!,/; ~ x i,l,k 


-B, 


i,l,k 


Vi,2,k yi,l,k 
^z,l,fc-f-l ^z,l,fc 
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(3.13) 


d 8 ~uj t i,jy 3 
dt 

di.i.JVa 
gi,l,JV 3 
ai.i.JVs 
e i,l,iV 3 
c i,l, 7 V 3 
<9^ u jj'i 

dt 

dij, 1 = 

bi,j,l ~ 
^J.l ~ 
e i,j,l = 

C *.j,l = 
= 

9S U i,j,N 3 

dt 

d i,j,N 3 


— di'i'XzS ~V?i-l,l,N3 + gi.l.JVs^ "^t.l.JVs-l 4 ^i,l,N 3 ^i,l,N 3 
+e iXN3 6^i +1 ,i,N 3 4 Ci,l,N 3 f>^? i,2,N 3 


A + 
A i, 1, 


iV 3 


%i,l,N 3 ~ x i- 1,1, N 3 

r+ 

■^i,l,N3 ^i,l,iV3— 1 


A- 


i.l.Afe 


41, 


n 3 


4 


B i,l,N 3 


r+ 

U i,l,N 3 


X i+1,1,N 3 “ Z»,1,7V 3 x i,l,N 3 ~ X i—1,1,N 3 Vi,2,N 3 ~ yi,l,N 3 z i,l,N 3 ~ z i,l,N 3 - 

~ A i,l,N 3 


X i+1,1,N 3 ~ x i,l,N 3 


B i,l,N 3 


Ui,2,N 3 ~ Vi,l,N 3 


di,j,lS U 4" U i,j — 1,1 4 U i,j,l 4 1<5 U l,j,l 

J r c i,j,l$~u 4 fi,j,\S~u j = 2 , * • • , N 2 — 1 


A + 

1 

x i,j, 1 — 1,^,1 

g a.i 

Vi,j, 1 1,1 


4 + 


4 


s .j,i 


afc,i 




C5.1 


*£i+l,j,l X i t j> 1 1,3,1 Ri+1,1 2/i, j, 1 yi,j,l yi,j — 1,1 z i,j,2 z i 

1 


-CTai 

4, 3, 2 ~ 


— diJ,Ar 3 ^ U i—l t j t N 3 4 b i,j,N 3 fi u i,j-l,N 3 4 g i,j,N 3 fi u i,j,N 3 -l 

^~&i,j,N 3 3 u i,j,N 3 4 ei,j,N 3 f> u i+l,j,N 3 4 c i,j,N 3 ^ u i,j+l,N 3 i j = 2, • • • , N2 ~ 1 

A + 

_ A iJ,Na 

x i,j,N 3 x i—l,j,N 3 


(3.14) 


,3,1 


(3.15) 
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\j,N 3 — 


Si,j,N 3 ~ 


*i,j,N 3 


R+ 

° i,j,N 3 


Vi,j,N 3 - yi,j~l,N 3 




z i,j,N 3 z i,j,N 3 — 1 


Ai,j,N 3 


A + 

A i,j,N 3 


+ 


B i,j,N 3 


x i+l,j,N 3 x i,j,N 3 x i,j,N 3 x i—l,j,N 3 Vi,j+\,N 3 yi,j,N 3 


B tj,N 3 


C tj,N 3 


e i,j,N 3 = 
c i,j,N 3 — 
d 8 ~Ui i N 2 ,l 


yi,j,N 3 yi,j—l,N 3 z i,j,N 3 z i,j,N 3 — 1 

~ A i,j,N 3 


x i+l,j,N 3 x i,j,N 3 


' B i,j,N 3 


yi,j+l,N 3 yi,j,N 3 


dt 


— di,N 2 ,l8 U i-l t N 2 ,l + bi,N 2 ,l8 U i,N 2 -l,l + &i,N 2 ,\8 U i,N 2 ,l 
+ e i,N 2 ,l8~U i+l,N 2 ,l + U,N 2 ,\8~^ i,N 2t 2 


H,N 2> 1 = 


>i,7V 2) l 


a;,iv 2 ,i = 


e i,N 2 ,l — 


H,N 2 ,1 


A + 

2,1 


x i,N 2t l x i~l,N 2 ,l 

B + 

° i ,JV a ,l 

2/i,Ar 2 ,l — yi,N 2 - 1,1 
A,7V 2 ,1 


4+ D+ 

A t,JV 2 ,l n i,N 2 ,l 


+ 


^t,W 2 ,l 


<9<5 "n 


i,N 2 ,k 


dt 


x i-\-l,N 2t l x i,N 2 ,l x i,N 2 ,l X i—1,N 2 ,1 yi,N 2 ,l yi,N 2 — 1,1 ^i,N 2 ,2 % i,N 2 , 

- A i,N 2 ,l 

X i+1,N 2 ,1 ~ x i,N 2 ,l 

~ (B i,N 2 , 1 
z i,N 2 ,2 — ^i,7V 2 ,l 

- di,N 2 ,k8^U i-l,N 2 ,k +^i,N 2 ,k8~U i,N 2 -l,k + gi,N 2 ,k8~U i,N 2 ,k-l 

+^i,N 2 ,k8~U i,N 2 ,k + ^i,N 2 ,k8^L i+l,N 2 ,k + U,N 2 ,k8~U i,N 2 ,k+l, k = 2, • • • , N$ 


H,N 2 ,k — 


H,N 2 ,k 


A + 

^-i'N'iyk 


x i,N 2 ,k x i—l,N 2 ,k 


B + 

U i,N 2 ,k 


Vi,N 2 ,k ~ yi,N 2 -l,k 


(3.16) 


(3.17) 
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g i,N 2 ,k = 
&i,N 2 ,k = 

&i,N 2 ,k — 

f i,N 2 ,k — 

d6T?i,N 2 ,N 3 

dt 

di.7V2.iV3 

ki,JV 2 ,7V 3 

gi.iVa.TVa 

ai,7V 2 ,7V 3 

e i,7V 2 ,iV3 


^? 7 Va.fe 


z i,N 2 ,k ~ z i,N 2 ,k - 1 


A, 7 V 2 ,fc 


4 + 

^i,N 2 ,k 


Bt,N 2 ,k 


^i+l,JV 2 ,fc - x i,N 2 ,k x i,N 2 ,k ~ ^i-l.TVa.fc 2 /i,A 7 2 ,A: " yi,JV 2 -l,fc 


+ ‘ 


c; 


i,N 2 ,k 


^i,N 2 ,k 


z i,N 2 ,k-\-l z i,N 2 ,k z i,N 2 ,k 1 


— AJ 


i, 7 V 2 ,fc 


x i+l,N 2 ,k x i,N 2j k 


-c. 


i,N 2 ,k 


z i,N 2 ,k+l ~ z i,N 2 ,k 


di, 7 V 2 , 7 V 3 <^i— 1 , 7 V 2 , 7 V 3 + bi y N 2 ,N 3 t>~U i,N 2 -l,N 3 + gi,JV 2 , 7 V 3 i,N 2 , 7 V 3 -1 
+a-i,N 2 ,N 3 $ u i,N 2 ,N 3 + e i, 7 V 2 , 7 V 3 <?> w i+l, 7 V 2 ,JV 3 


( 3 . 18 ) 


A + 

^i.T^.TVs 


x i,N 2 ,N 3 ~ X i— 1 ,N 2 ,N 3 


n 2 ,n 3 


yi,N 2 ,N 3 ~ yi,N 2 -l,N 3 


Ct,N 2 ,N 3 


z i,N 2 ,N 3 ~ z i,N 2 ,N 3 — l 

a: 


\n 2 ,n 3 


A + 

A i,N 2 ,N 3 


x i+l, 7 V 2 , 7 V 3 ~Xi,N 2 ,N 3 Xi,N 2 ,N 3 ~ ^i-l,JV 2 ,N 3 


n i,N 2 ,N 3 


^i^N 2 ,N 3 


yi,N 2 ,N 3 ~ yi. 7 V 2 -l.JV 3 z i,N 2 ,N 3 ~ z i,N 2 ,N 3 -l 


-a: 


i| 7 V 2 , 7 V 3 


®i+l, 7 V 2 ,JV 3 — x i,N 2 ,N 3 
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Wall Boundary Conditions 



Equation 3.13 


Equation 3.15 


Equation 3.18 


3.2.2 Inflow Boundary Conditions for the System Matrix 


The inflow boundary conditions are implemented as modifications to the system matrix following the 
method developed in section 2.3.2. For this paper, the inflow boundary condition is assumed to be supersonic; 
therefore, there are five physical boundary conditions and zero numerical boundary conditions [10]. From 
equation (2.49) , the matrices for the supersonic inflow boundary condition become: 


Li 

L 2 


dB bc 

d~u 

0 


(3.19) 


and 


_ / — 1 D 

dt ~ L i B bc 

dB bc 




B bc 


with 


B bc = 


P 

m\ 

m 2 

m 3 

£ 


Equation (3.20) becomes the following when it is linearized: 

-i -l 


dSu 

dt 

If there are no upstream inputs 6B bc becomes: 


dB bc 


du 


SB bc 


_ dB bc 

0 i$bc — r\ V ® ^ 

« a u 

and the inflow boundary condition is represented by the following equation: 

dS U 


dt 


— ^ U lj'k 


with j = 1 , • • • , N 2 and k — 1 , • • • , A 3 . 


(3.20) 


(3.21) 


(3.22) 

(3.23) 

(3.24) 



Inflow Boundary Conditions 


These boundary conditions are included as modifications to the system matrix at the first grid point in the 
x-direction. 
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3.2.3 Outflow Boundary Conditions for the System Matrix 

The outflow boundary conditions are also implemented as modifications to the system matrix following 
the method developed in section 2.3.2. In general, the outflow boundary condition is subsonic; therefore, 
there are four numerical boundary conditions and one physical boundary condition [10]. The equations for 
the subsonic outflow boundary condition are as follows: 


dS u jv - ! 1 1 j c ► c > ¥ ^ 

02 ’ ’ = “^,1,10 u jvi-i,i,i + a n u i,io u Nu i t i -f cjVi.1,1 6 u Nl>2> \ + f/v 1 ( i ( i<5 u Nu1}2 ( 3 . 25 a) 


dSll N u l,k , c — >• , > r __> 

~ 7^ — O-Nul.kO u 7Vi-l,l,fc + gAT 1 ( l,fcC> U N u l,k-1 + SLN lt l,kO U N u l,k 

+CN u l,kS ~UN lt 2,k + fATi.l.fc^ATx.l.fc+l, fc = 2, • ■ • , 7V 3 — 1 


(3.25b) 


d5 U t \ t N 3 J C - Tr y c I r , 

Q t ' ’ = &Nx,\,N 3 0 u JVi — l.l.ATa + gN u l,N 3 fi u N lt l,N 3 -l + ^N lt l,N 3 8 U N u l,N 3 + C Nl ,l,N 3 ^ N lt 2,N 3 

(3.25c) 

d8~u c _ x y r 

^ = d N u j,i8 u + biVi,i,i<5 u N t j- 1,1 + a^! J,1<5 u ^,^1 

-few! j+i,i + fjvij.i 8~u j = 2, • • • , N 2 — 1 (3.25d) 


d8~U N u j,k 

dt 


d N lf j,k8ll ^-lj.k + b Nl ,j,k8~U N u j-l,k + SN u j,k^N u j,k-l 
+ a -N 1 ,j,k8~U Ni,j,k + C Nl ,j,k8ll Ni,j+l,k + ^N lt j,kSll N lt j,k+1 

j = 2, - • • , AT 2 - 1, *r = 2, - - - , JV 3 — 1 


(3.25e) 


d8~u 


N u j,N 3 


Q t — ^N u j,N 3 8 u Nl - l,j,N 3 + ^>Ni,j,N 3 8 U N lt j-1,N 3 + gNi ,j,N 3 81t N lt j,N 3 -l 

+ a N l) j,N 3 8'u N u j,N 3 + c N lt j,N 3 Sll N lt j+1,N 3 , j = 2, • • ■ , N 2 - 1 (3.25f) 


d8~u Nx,n 2 , l 
dt 


— djV 1 ( jV 2 ,l<5 U ATi-l.JVa.l + bjv ll JV 2 ,l 6 '^jV 1 ) ./V 2 -l,l + ajv^jVj.l^ "^Ar^iVa.l +fAr 1 ,7V 2 ) l<5u > 7V 1 ,iV 2l 2 

■ (3.25g) 


ddll N x t N 2 ,k 

dt 


^N u N 2 ,k81l jVi_l.JVa.fc + biv ll jv a ,fe5't?jv 1 ,jv a _i f jb + gN!,N 2 ,k8^U N u N 2 ,k-l 

+ a 7v 1 ,N 2 ,k8~U Ni,N 2 ,k + fNi,7V 2 ,fc^W > Ar 1 ,iV 2( fc+l, A? = 2, • • • , AT 3 — 1 (3.25h) 


951? n x ,n 2 ,n 3 
dt 


^N U N 2 ,N 3 ^ JVi-l.JVa.iV 3 + bjVuTVa.iVa^ l^JVi.iVz-l.TVa + ZN 1 ,N 2 ,N 3 8'u N u N 2 ,N 3 -1 

+ a -N l ,N 2 ,N 3 8'u n u n 2 ,n 3 (3.25i) 
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Outflow Boundary Conditions 


• 

# 

• 

Equation 3.25a 

Equation 3.25b 

Equation 3.25c 

€> 


4 

Equation 3.25d 

Equation 3.25e 

Equation 3.25f 

• 

* 

• 

Equation 3.25g 

Equation 3.25h 

Equation 3.25i 


These boundary conditions are included as modifications to the system matrix at the last grid point in the 
x-direction. They will be discussed at length in section 3.3 where the coefficients of the equations are defined, 
and the input matrix, B, is derived. 

3.3 Input Matrix for Downstream Mach Number 

When the flow at the compressor face is subsonic, there are four numerical boundary conditions and 
one physical boundary condition [10]. The numerical boundary conditions are associated with the positive 
characteristics, and the physical boundary condition is associated with the negative characteristic. The 
implementation of the downstream Mach number as a boundary condition input is derived below. 

To begin with, take the inverse of (2.25) , and then partition the matrix following the procedure from 
section 2.3.2. 
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Therefore, for the downstream Mach number boundary condition, L\ is defined as: 


Li = 


dB bc 

dli 


(3.32) 


and L 2 is defined as: 


L 2 = 


— 1 \N 


(^r 1 ) 


0 


(3.33) 


From equation (2.48) , the compressor face boundary condition may be written as follows. 


d~u 


Ni ,j,k . r _! 


dt 


+ L x L 2 An u j t k 


dli 


Nij,k 


dx 


+ ,j,k 


dli 


Ni ,j,k 


dy 


+ C N lt j,k 


dli 


Ni,j,k 


dz 



0 

1 1 

— B bc 


(3.34) 


For this boundary condition to be implemented in the small perturbation model, it must be linearized as 
shown below. 


d6~u N u j,k 
dt 


I r — 1 r f a U N\,j,k . 7 -, 

+ L 1 L 2 I ,j,k f - B Nl J tk 


dSli Nu j, k , ^ d8~u Nl j k 4 r 

+ ^N u j,k " I = B 


dy 


dz 


-1 

1 


0 

-8B bc 


8B bc is calculated as: 


and can be rewritten as: 


co dB bcc-r+ , dB bc C1t , 
OBbc — o— > ^ U "b oar oMi n p U t 

d U dMinput 


SBbc — 


dM dM dM dM dM 
dp drrii dm2 dm3 de 


8 U 8 Ad input 


dS~\i 

Now substituting equation (3.37) into equation (3.35) and replacing with equation (2.31) yields: 

dx 


(3.35) 

(3.36) 

(3.37) 


d8~u Nj ,j,k 
dt 


= —L7 l LoA 


1 ,j,k 


fi~UNi,j,k ~ 8 ~Un 1 — l,j,fc 
x N lt j,k ~ x Nx~l,j,k 


— 1 r d d8~u Nltjt k ! r d81i Nltjtk 

~ L 1 B 2 L>Nx,j,k' 


—L\ l L 2 B Nltj>k - 


dy 


—L 


-1 


+Ti 


-1 


04x5 

dM dM_ dM_ dM_ dM 
dp dm\ dm2 dm3 de 
-1 T 


dz 


8~u 


Ni,j,k 


(3.38) 


0 0 0 0 1 


8 Mr 


input 
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this can be rewritten as the following: 


d8!?Nu3,k 

dt 


= £r%(- 


c—’ , ^Ni,i,k 

0 u N, ,j,k H 


x N lt j,k ~ XNi-l,j,k 

, °N U j,k CTT> 

H 0 W JViJ-l.fc 

yN u j,k - yN u j-i,k 

B 


x Ni,j,k X N X — l,j,k 

B+ 

~Nujjk 


■81lNi~l,j,k 


B 


N u j,k 


yN lt j,k - yN u j-l,k VNiJ+l,k - VNi ,j,k 


n+ 

N lt j,k c — ♦ , ° N u j,k c — > 

0 u N u j+l,k + — — 0 U N lt j,k- 1 


yjvu'+i.fe - yATi.j.fc 


- z N\,j,k—l 


^ Ni,j,k 


(3.39) 


C 


Wi.j.fc 


z Ni,j,k z Ni,j,k—l z Ni,j,k + 1 z Ni,j,k 


8~U 


c 


Ni ,j,k 


N u j,k 


z Ni,j,k+l z Ni,j,k 


■8~U 


Ni ,j,k-\- 1 


-^r 1 


04x5 

dB bc 

L <91? 


<5 u JV t ,j,fc + 1/! 


-l 


-|T 


0 0 0 0 1 


8M, 


input 


The terms for the input matrix, B, are obtained from the coefficient on 8M input . Equations (3.40) through 
(3.48) are used to modify the system matrix, (3.10) . 


d8~u jVi ,i,i 


dt 

d/^,1,1 = 


U Ni- 1,1,1 + &Ni, l,l8~U Ni, 1,1 + CATi,l,l^l?7Vi,2,l + fjv lf 1,1^1?/^!, 1,2 (3.40) 

■4iVi,l,l 


ajVi.i.i = L x l L,2 

CATi.1,1 = —Li l L,2 

f/Vi ,i,i = —L 1 1 L2 


x Ni,1,1 — X N\— 1,1,1 

Ani,1,1 


i \- 


x Ni,l,l - a?JV! -1,1,1 VN lt 2,1 - yiVi.I.l ^,1,2 - z Ni ,1,1 


— L 


-l 


04x5 

dB bc 


L <9' 


B 


Nj.,1,1 


yN u 2,l ~ VN U 1,1 

QnTi,1,1 

2 iVi,l,2 “ *JVi,l,l 


d8~u ,i,fc 

dt 


&N lt l,k&~rf Ni-l,l t k +SN u l,k8~U Nul,k-1 +ajv 1 ,l,fc^l?JV 1 ,l,fc 
+CNi,l,fc«'u > JV ll 2,fc + fjV lf l,fc^l? JVi,l,fc+l j fc = 2, ■ • • , A^3 - 1 


(3.41) 


div 1( i,fc = Tj 1 L>2 


A-N u l,k 


x N lt l,k ~ x Ni- 1,1, fc 


gN u l,k = X L2 


r » + 


z N\,l,k ~ z N u l,k-l 
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a/Vi,i,fc = L x l L*2 — 


An u \, 


1 ,k 


+ 


B 


Ni,l,k 


x N u l,k ~ XNi-1,1 ,k yN u 2,k ~ yN lt l,k 


^N u l,k 


+ 


c 


N it l ,k 


z Ni,l,k ~ z Ni,l,k-l z N lt l,k+l ~ z N u l,k 


-LI 1 


C>4x5 

dBbc 

L dli J 


c N u l,k 

^N u l,k 


= -l^l 2 


= -LfLi 


B 


N l ,l,k 


yN lt 2,k ~ VNi,l,k J 

L'Nul.k 

z N lt l,k+l ~ z N u \,k 


d81l N 1 , 1 , jy 3 

dt 



= ir 1 ^ 

gJVx.l.JVa 

CS 
7 r-4 

II 

ajv lt i,7v 3 

= L^L 2 


d Ni,\,N 3 & u Ni-1,1,N 3 + SNi, 1,N 3 ^ N lt l,N 3 ~l 
+^N lt l,N 3 6 1? N lt l,N 3 + c N 1 ,l,N 3 SllN 1 ,2,N 3 

ANi,1,N 3 


x N u 1,N 3 — ^N!-l,l,N 3 

r+ 

^N u 1,N 3 

z N u l,N 3 — Z N\ , 1 , N 3 — 1 

Ani,1,N 3 


+ 


'Ni,1,N 3 


{ x Ni, 1,N 3 - X N 1 -1,1,N 3 yN lt 2,N 3 -yN u l,N 3 

Cn u 1,N 3 


CNi,1,N 3 = —L l 1 L,2 


z Ni,1,N 3 — Z N X , 1,N 3 -1 

B\ 


-LI 1 


04x5 

dBbc 

L dli J 


'Ni,1,N 3 


yN u 2,N 3 —yNi,l,N a 


d6~u 


, j, 1 


dt 


d N lt j,l = 




d N u j,iS u + h Nltj>l 6!? Nuj _ ltl + a-N u j,iS~u N^j.i 

+C Nl ,j,lSll Nij+1,1 + fN lt j,l$ll N 1} j,2i j = 2, • • • , jV 2 - 1 


li 1 l 2 


A 




— L l 1 L 2 


l^l 2 i - 


x Ni,j,l x Ni — l,j,l 

R+ 

U Ny,j,l 

yN u j,i -yN u j- i,i 

A N u j,l 


%,J,1 


_ yJVi.j.l - 2/iVi.i-l,! 


B 


N u j, 1 


+ 


c 


Wi .j'.l 


2/M J+1,1 ~ yNi,j,l z N u j,2 - z N lt j,l 


LI 1 


04x5 

dBb c 

L < 9 "u J 


(3.42) 


(3.43) 
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c N u j, 


1 = -L^L 2 




VNuj+1,1 -VNuj, 1 

' ^,1 ^ 
z N 1 ,j,2 - z Nu j, i / 


d8~u Ni,j,k j ► ^ ^ 

’ ’ = d N u j,k u Nv-lJ.k + b N lt j,k u N u j-l,k + g N u j,k u Nuj.k - 1 + &Nij,k V N u j,k (3.44) 

+ c iVi,i,fc^iV l ,i+i,fc + iVLi.fc+1, 3 = 2, • • • , N 2 - 1, k = 2, • • • , jV 3 - 1 


d Nx,j,k 

— Tr l r o 1 

f ,j,k 

-L /2 | 

K X Nx ,j,k — %Nx— l,j, k 

biVi,j,fc 

— rr l r,o 1 

( B N u j,k 


[^yN u j,k-yN lt j-i,k 

SN u j,k 

— rr 1 1 

( C Ni,j,k 

-^1 -^2 

\^ Z Ni ,j,k ~ z Ni,j,k-l ^ 

a Ni,j,k 

= £r%| 

( A-N u j,k 


a NlJ , fc = L^LJ + 

y x Ni,j,k x Ni — l,j,k VNuj.k VNi,j — l,k VNi,j+l,k VNi,j,k 

C N } ,j,k 


B 


N u j,k 


I 

z N\,j,k z Ni,j,k—l z Ni,j,k+l z Ni,j,k 


L\ 


04x5 

dB bc 

d~u J 


C Ni,j,k — L y L 2 


B 


Nu3,k 


\ 

\,j,k — —L\ 1 L 2 f 


yN u j+i,k - yN u j,k 


f AT, 


dt 


C 


Ni,j,k 


d N u j,N 3 

II 

hT 

to 

b N u j,N 3 

= Lf'Lj 

SN u j,N 3 

= Li'Ll 

a N L ,j,N 3 

= L^L 2 


\y Z Ni,j,k+l z Ni,j,k J 

d N lt j,N 3 81? Ni-l,j,N 3 + b^ J,_/V 3 8~U Ny + &N u j,N 3 8ll Ni,j,N 3 -l 
+a. Nl> j t N 3 8' ~UN u j,N 3 + c N u j,N 3 81l Nl ,j+1 )N3 , j = 2, • • • ,N 2 - 1 

A N lt j,N 3 

, x N u j,N 3 ~X Nl -l,j,N 3 

f R+ 

~Ni ,j,N 3 


yNuj.Na ~ yN lt j-l,N 3 

n+ 

U jVi,j,N 3 

z Ni,j,N 3 z Ni,j,N 3 — l 

< J 

A Nx,j,N 3 


B+ 

n Nx,j,N 3 


x N l} j,N 3 ~ x Ni~l,j,N 3 yN u j,N 3 ~ VNx,j- 1,N 3 


+ 


B N u j,N 3 


r+ 


yNi,j+l,N 3 yNi ,j,N 3 Z N\ ,j,N 3 Z N\ ,j,N 3 — 1 


-L 


-l 


04x5 

dB bc 
L dli J 


(3.45) 
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c N lt j,N a 


-l^l 2 ( 


VNuj+l.Ns -VN lt j,N a 


dSll N lt N< 2 ,l 
dt 


dN u N 2 ,lfi u N x - 1,N 2 ,1 + ^Ni,N 2 , 1<S U N x ,N 2 ~\,\ 

+^N x ,N 2 ,lS ^JVj.JVz.I + fWi.TVa.l^JVi.A^ 


- li x l 2 


bjvi,jv 2 , i = L^L 2 


ANi,N 2 ,1 


a JVi,JV 2 ,l — L x l L 2 [ — 


x N lt N 2 , 1 - x Ni~l,N 2i l 

B + 

n N u N 2 ,l 

yN u N 2 ,l y N x t N 2 — 1,1 

An x ,n 2 ,i 


R+ 


x N\,N 2 ,1 ~ %Ni- 1 ,N 2 ,1 yN u N 2 ,l ~ yNi,N 2 -l,l 


+■ 


C 


JVi,JV 2 ,l 


z N x ,N 2 ,2 ~ ZN u N 2 ,1 


LT l 


04 x 5 

dB bc 

L dli J 


^Ni,n 2 ,i — —L x x L 2 


C 


n u n 2 ,i 


z N u N 2 ,2 ~ z N lt N 2 ,l 


dSll N u N 2 ,k 

dt 


djVi,JV 2 ,fc 

<N 

u* 

II 

b7V1.Ar2.fc 

II 

Jpl 

h-* 

IO 

gN1.JV2.fc 

= ir% 

a JVi,JV 2 ,fc 

= 


d;Vi,W 2) fc<5 U iVi-l ,N 2 ,k + b N lt N 2 ,k 61 ? Nx,N 2 -l,k + gN u N 2 ,k 8 ll NuN 2 ,k-l 
d-^N x ,N 2 ,kSll N x ,N 2 ,k + ^N lt N 2 ,k^ ~U N x ,N 2 ,k+l, k = 2 , ■ • • , JV3 — 1 

ANt,N 2 ,k 


,N 2 ,k ~ l,JV 2 ,fc 


R+ 

n N u N 2 ,k 


yNi,N 2 ,k ~ yNi,N 2 -l,k 


r+ 

0 Ni,N 2i k 


z N u N 2 ,k ~ ZNi,N 2 , fc-1 
A/Vi.JVa.fc 


R+ 

^ATi.JVz.fc 


x N u N 2 ,k ~ ffJVi-l.Afe.fc yN lt N 2 ,k ~ yN u N 2 -l,k 


f~*+ 

<w ''7Vi,7V 2> fc 


+ 


c 


N u N 2 ,k 


-1 


z N u N 2 ,k ~ ^AT lf JVa,fc-l ,JV 2 ,fc+l ~ z N y ,N 2 ,k 


04 x 5 

dBbc 

L dl? J 


fjVi,7V 2 ,fc = — 1 L2 


C 


N lt N 2 ,k 


z Ni,N 2 ,k+l ~ z M x ,N 2 ,k 


(3.46) 


(3.47) 
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dS U N t ,N2,N3 

dt 


d-Ni,N 2 ,N 3 

H 
1 H 

II 

bjV1.JV2.JV3 

= Lph 

Sn u n 2 ,n 3 

= LpLi 

*N lt N 2 ,N 3 

= LpL 2 


&N lt N 2 ,N 3 &~rf ATi-l.TVa.ATg + bN u N 2 ,N 3 ^Ni,N 2 -l,N 3 
+%Ni,N 2 ,N 3 6 U N u N 2 ,N 3 - 1 + a iVi ,N 2t N 3 ^^ N lt N 2l N 3 

Ani,N 2 ,N 3 

X N U N 2 ,N 3 ~ x N 1 -1,N 2 ,N 3 
R+ 

n Ni,N 2 ,N 3 


yN u N 2 ,N 3 ~ VNuNi-l.Na 
' sv+ 

^Ni,N 2 ,N 3 

Z N X ,N 2 ,N 3 ~ ZNi,N 2 ,N 3 -1 

\ 

An u n 2 ,n 3 


B + 

°N u N 2t N 3 


x N\,N 2 ,N 3 — x Ni~1,N 2 ,N 3 VNi,N 2 ,N 3 ~yN lf N 2 -l,N 3 

-l; 


Cn u n 2 ,n 3 


z Ni,N 2 ,N 3 - ZN u N 2 ,N 3 -1 


--1 

J 1 


04 x 5 

dBhc 


L d~iL \ 


(3.48) 


3.4 Output Matrix 

The outputs for the linearized inlet model are static pressure, total pressure, and Mach number. For an 
input located at the compressor face, these outputs are located downstream of the normal shock. Outputs 
upstream of the normal shock remain unaffected by a downstream disturbance, therefore they do not need 
to be considered. 


3.4.1 Static Pressure 


The components of the output matrix, C, for the static pressure response are obtained by linearizing the 
static pressure equation: 


V = (7 - 1) ( e “ 


mf + 1712 + m 3 

2 p 

dp dp dp ,9p dp 

op = -^-6p+- — 8mi + - — 8m 2 + - — < 5 m 3 -f — 8e 

dp dmi dm 2 dm% de 


8p — 


(3.49) 


5 (7 - 1) ( mi+ 7 + m3 ) - (7 ■ - 1) (^) ^ - (7 - 1) (^) *r» 

- (t - 1) + (7 - 1) 

Each coefficient of the small perturbation terms is an entry in the output matrix. 


3.4.2 Total Pressure 

The output equation for the total pressure is determined in the same manner as for the static pressure. 
The components of the output matrix for the total pressure response are obtained from the linearization of 
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the total pressure equation: 


Pt = p ( 1 -f 


7-1 ' |f ~ 1 


(3.50) 


The pressure and Mach number terms are replaced with expressions in terms of p, mi, m 2 , m3, and s, 

2£7 P + (mf + m| + m§) (1 - 7) \ ' > '“ 1 


(7 - 1) (2s7p + (mf 4 - m| + m§) (1 - 7)) 


Pt = 


7 (2s p — mf — m| — m|) 


27/? 


(3.51) 


(3.52) 


Then the total pressure differential is calculated as shown below. 

c dpt r . dpt dpt dpt dpt _ 

opt = -x— <5/?+ a — om\ + ^ — v fn 2 + -X — Sms -f -z-Se 
dp ami dm2 dm$ ds 

Each one of the coefficients in the small perturbation terms of the above equation is an entry in the output 
matrix for the total pressure. 

dpt = ^ ( 2 ^ (T ■ ~ 2) - (7 2 - 47 + 1)) (- ^/ £ l * 2,~ 7) ) 7 ~ 1 

dp 27 (2 s — p~v 2 ) 

1 


dp t 
dm 1 


dpt 

dm2 


dpt 

dm 3 


dpt 

ds 


U (^ - 37 + 1) - 2e 7 (7 ~ 2)) - p ^> 2) 7) ) 7 1 

7 (2s — pi? 2 ) 


V ( P V 2 (T» - 87 + 1) - tor (7 - 2)) (- 2 ^|^ ! 7) ) 7 1 

7(2£-pV 2 ) 

™(pV 2 ( 7 2 -3 7 + l)-2 e7 ( 7 -2)) 

7 (2s — pT* 2 ) 

(2 £ 7 (t - 1) - pV 2 (T 2 - 7 + 1)) ( 2£7 7 | 2 ^ ” 2 p ^~ 7) ) 7 ~ 1 

7 (2s — pi? 2 ) 


(3.53) 


3.4.3 Mach Number 

Likewise, the output equation for the Mach number is determined in the same manner. The components 
of the output matrix for the Mach number response are obtained from the linearization of the Mach number 
equation, 

m= V^TZ±Z (3.54) 

The velocity and speed of sound terms are replaced with expressions in terms of p, mi, m 2 , m3, and s (as 
shown in section 3.3). Again, the total differential of the Mach number is calculated, 


,,, dM. dM . dM dM c dM . 

oM = — op + - — dmi + - — om 2 + 7; — < 5 m 3 + ^r-< 5 s 
ap ami am 2 dm3 ds 


(3.55) 
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and each coefficient on the small perturbation terms is an entry in the output matrix for the Mach number. 


m 

dp 


dM 

dm\ 


y/2 y/m\ + 4 - (e — pT* 2 ) 


'y ( 2 e (7 — 1 ) — 7 p~v 2 ) 4- p~v 2 


y(2£-p~v* 2 ) (1 — 7 ) |p| 


~ ^ ~ 7ply>2 ) + p ~^ 

7 ( 2 e — pT* 2 ) (7 — 1) yj m\ 4 - m 2 + m\ 


sgn(p) 


dM 

dm 2 


dM 

dm 3 


\/2m 2 ‘ 


1 7 (2e (7 — 1) — 7 pU >2 ) 4- p~v 


sgn (p) 


7 ( 2 e — pH* 2 ) (7 — 1 ) yj m 2 4 - m| 4 - ra§ 


y /o m ^ ^ 7 ( 2 £ (7 ~ 1 ) ~ yp^v 2 ) + P^ 

7 ( 2 e - pV 2 ) (7 - 1 ) yj m\ + m| 4 - m§ 


sgn (p) 


(3.56) 


/— 5 5 5 / 7 (2e (7 - 1) - yp~v 2 ) 4- pV 2 

dM V2y/mi+mi + miJ — l sgn (p) 

< 9 e 7 (2e — pl^ 2 ) 2 (1 — 7) 


4 Model Reduction 

The models developed from the linearized split flux method are usually too large to be used effectively 
in the design of a control system; therefore, model reduction is necessary so that the linear model can be 
transformed into one that is of manageable size. In addition, a reduced order model (ROM) is needed when a 
transfer function of the system is desired, because the calculation of the transfer function from the full order 
model (FOM) becomes a numerically ill-conditioned problem. In general, the linear reduced order model is 
expected to perform like the linear full order model. There are many different methods that can be used for 
model reduction; those considered here take advantage of the state space format of the linear model. 

One of the most common model reduction methods is singular perturbation [13]. This method requires 
a linear model that can be partitioned into slow and fast subsystems so that a reduced order model can be 
obtained by neglecting the fast subsystem. The ability to partition the model into subsystems indicates that 
the model possesses a two-time scale property; that is, there is a large gap in the spread of the eigenvalues 
of the linear model. There is little contribution to the dynamics of the system from the fast eigenvalues; 
therefore they can be discarded and a reduced order model obtained. But as was shown in [14] this method 
does not provide the smallest reduced order model that can be achieved, since in general the linear systems 
do not possess a two-time scale property. 

Another popular method is balancing [15] [16] [17]. This method requires a linear model that can be 
partitioned into a strongly controllable/observable subsystem and a least controllable/observable subsys- 
tem. A reduced order model is then be obtained by discarding the least controllable/observable subsystem. 
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However, the ability to partition the system into these subsystems indicates that the model is minimal; 
that is, it is both controllable and observable. In general, linear models developed from the linearized CFD 
method are not minimal. Therefore this method does not work well because of the many uncontrollable and 
unobservable states in the linear model. 

Two model reduction methods that are successfully applied to the linear model are the Schur and square 
root methods [18]. Both of these methods take into account the controllability and observability of the Unear 
system, but they do not depend on the Unear model being minimal. Although the Schur method is more 
numerically robust than the square root method, the square root method is preferred over the Schur method 
because it is computationally less expensive. In addition, since the reduced order model is balanced from the 
square root method, smaller reduced order models can be extracted from the original reduced order model 
without any more computations. The same is not true for the Schur method. The square root method is 
described in the next section and will be applied to the linear models developed later in the paper. 

4.1 Square Root Model Reduction 

The square root method of model reduction detailed in [18] is summarized here for the readers convenience. 
This model reduction method reduces the full order model by a balancing transformation that requires only 
the first k (size of reduced order model) columns of the balancing transformation matrix. The reduced order 
model (Ak, Bk, Ck, Dk) for the state space system (A, B, C, D) is calculated from: 

A-k S ' R,big 

B k = S T LMg B 
Ck = CS ft t big 

D k - D 

Computation of the transformation matrices S ^ big and S^big is shown below. 

To begin with, the controllability and observability Grammians, P and Q, are calculated from the fol- 
lowing Liapunov equations. 

A P + PA t = -BB t 

A t Q + QA = -C r C (4.2) 

Once the Liapunov equations are solved for the Grammians, the square roots of P and Q ) 

P = L c L t c 

Q = L 0 L t (4.3) 

are calculated from the following: 

[U P ,S P ,V P ] = SVD(P) 

L c = U p diag 

(4.4) 

iu g ,s qt v q ] = SVD(Q) 

L 0 — Uqdiag 

A singular value decomposition ( SVD ) of L^L C produces the matrices needed to compute S£ big and Sji tbig . 

\U^V\ = SVD{L t 0 L c ) (4.5) 




(4.1) 
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^ is defined as: 


El = diag(ai,(T 2 , • * *, <An) (4.6) 

where m is the number of nonzero Hankel singular values represented by a. The a are defined as the square 
roots of the eigenvalues of the PQ matrix. 


<7 = (A (PQ))* 


(4.7) 


Prom Sfcai 1S defined as: 

E = diag(cri,cr 2 , • • •, cr k ) (4.8) 

Note that the a are in descending order according to their magnitude. Now S^big an d Sr^ can be calculated 
from the following equations, 


$L,big 
Sr, big 


= L 0 U 


= L C V 


r 

E bed 
0 

Yhbal 

0 


(4.9) 


and the reduced order model can be determined using equation (4.1). While this method works well for 
the linearized one dimensional CFD models, it must be modified before it can be used for the linearized 
multidimensional CFD models; because due to their large size and sparse structure, the multidimensional 
CFD models are numerically ill-conditioned. These modifications are described in the next section. 


4.2 Modified Square Root Model Reduction 

Unfortunately for the linearized multidimensional CFD models, the calculation of the Grammians that 
are needed for the square root model reduction is computationally prohibitive because of the large number 
of equations involved. However, low-rank approximate solutions to these large scale Liapunov equations are 
possible. The approximate solution for the large scale Liapunov equation is calculated using the Arnoldi 
method from [19]. In addition, further numerical savings are achieved if the SVD calculations are performed 
using the partial SVD method from [20]. The modifications that are needed to the square root model 
reduction method when using these approximations are detailed in this section. 

Since the controllability and observability Grammians are approximations, the square root model reduc- 
tion method has to be modified [21]. The approximate controllability and observability Grammians are equal 
to: 


approx 


Q approx 


= V„x p v . P T 

= V q X,V? 


(4.10) 


where the matrices V p , X p , V q , and X q are calculated using the Arnoldi approximation method from [19]. 
Then the matrices L c and L 0 from (4.3) are calculated using the following equations: 


Lc — VpRp 

L 0 = V q Rq 


where Ftp is the square root of X p , and R q is the square root of X q as shown below. 


[Ux p ,S Xp ,V Xp ] 

Rp 


SVD (X p ) 

Ux 9 diag {^Jdiag (S Xp )^j 


(4.11) 
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(4.12) 


[Ux„Sx„V Xq ] = SVD(Xg) 

Rq = V Xq diag (^Jdiag (S x ,)j 

The matrices needed to compute S^ big and Sr^q for equation (4.9) are then calculated from: 

{V,'Ei,V) = SVD(lZL c ) (4.13) 


where 1S defined (4.6) by and J^bai defined by (4.8). Finally Shjbig an d SrMq are determined from 
the following equations, 


SL,big 


^R, big 




L C V 



(4.14) 


and the reduced order model is computed using equation (4.1). For this work, the singular value decomposi- 
tions in the above equations were calculated using the subspace approximation method from [20] . Note that 
in the modified square root model reduction, the approximate Grammians are never calculated; this results 
in a significant savings in computation time and in storage requirements. The errors associated with the 
model reduction and linearization processes are described in the next section. 


5 Uncertainty 

There is a certain amount of error associated with the linear models developed using the linearized split 
flux method; in particular, there is error due to the linearization process and the model reduction process. 
Tkese errors can be characterized as modeling uncertainties which are represented by error bounds and are 
derived below. 


5.1 Linearization Error 

The major source of error incurred in developing linear control models lies in the linearization process. 
As long as all the expected perturbations from steady-state are “small”, the error bounds will be small. 
Here, small depends upon which variables are being considered and how sensitive the system dynamics are 
to steady-state changes. Clearly this must be determined by trial and error for each system considered. 
Normalizing the data allows all the variables to be weighted similarly. While many approaches can be 
considered for determining error bounds, the one described below uses infinity norm bounds and is probably 
the most useful for the control system designer. In addition, it is also a convenient form for the numerical 
experiments considered. It should be noted that this error bound holds exactly only for linear systems; as 
such, it only serves as an approximation of the true bound. Its accuracy will certainly be degraded for larger 
perturbations. 

The modeling error due to the linearization process is based on the fact that the infinity norm of the 
transfer function error is less than or equal to the one norm of the unit impulse response error (B. Veillette, 
University of Akron, Akron, OH., private communication). 

]|Gt rue CM - G(ju/)lloo ^ ligtrue ( £ ) “ g( £ )ill C 5 ’ 1 ) 

Gtme (jM represents the transfer function of an exact linear model, if it were available, and G (jco) represents 
the transfer function from the split flux linearization process. The infinity norm of the transfer function error 
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is defined as: 

||Grf rlxe (jo>) G (j^’)|| CX ) — ^rnax {G^ rue (jtu) G (j'cu)} 

and the one norm of the impulse response error is defined as: 


J rOO 

|S true (p) — S Wl dt 

0 

fo° |g true (t) — g {t)\dt can be approximated by: 


(5.2) 


(5.3) 


N 

|g true (t) ~ g (i)J dt ~ \Strue “ g”| AT (5.4) 

n — 0 

Therefore, the error bound on the infinity norm of the transfer function can be represented by the following: 

N 

||G( rue (ju) - G OHIL < X] |g?,„e - g"l AT (5.5) 

n=0 

The notation can be simplified by expressing G trU e (jv) — G (jio) as E (ju>) and g^ rue — g n as e n . Therefore, 
the linearization error bound is represented as: 

N 

l|E(ju/)IL<5>”1AT (5.6) 

For multi-input/multi-output systems (MIMO), the transfer functions and impulse responses are represented 
in matrix form, and the linearization error bound can be determined from the following: 



|E(* 


lien Will 

lleiaWlli * 

• IhipWII, 

Hear Will 

IleaaWlli ‘ 

||®2p (p) || 1 

||e r i (t ) || 1 

lleraWlli • 

||e r p(t)||i 

is defined as: 




(5.7) 


a max = (A max (matrix 7 * matrix)) 


(5.8) 


For a MIMO system with p number of inputs and r number of outputs, the transfer function error matrix is 
written as: 

En (ju) E12 (. ju) ■ ■ • Ei p (jca) 

E21 (jw) E22 (ja>) ■ ■ ■ E2 p (jui) 


E = 

L Erl 0‘w) E r2 (juS) 

and the impulse error matrix is written as: 


E rp (ja>) J 


(5.9) 


en (t) 

ei2 (t) ■ • 

e ip (t) 

e2i (t) 

e22 (t) • ■ 

• • e 2p W 

Grl (t) 

e r2 (t) • ■ 

e rp (p) 


(5.10) 
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In the preceding definition, the one norm of each element is: 


rOO 

/ |e| dt — 

Jo 


and the approximation for (5.11) is: 


f 


lei dt 


lien (t)\\i 

IteWlli ■■ 

'• U e lp(*)Hl 


Ileal WII, 

ll®22 (l)lli '■ 

|| e 2p(t)|| 1 

(5.11) 

. KiWII, 

IteWlli 

1 

T“t 

0 


N 

N 

N 

1 


E left | AT E Mil AT 

n=0 n=0 

E tel at e |e&| at 

n — 0 n—0 


E tel AT E tel AT 

L n—0 n=0 


E tel AT 

n=0 

E K|AT 

n=0 


E K|A T 

n=0 


(5.12) 


Therefore the infinity norm of the error transfer function matrix is bounded by: 


il E C?w)[loo 



N 



N 

N 





E 

e ii 

AT 

E tel at ■ 

• E 

| e ip 

AT 



n — 0 



71=0 

71=0 




N 



N 

N 





E 

e 21 

AT 

E tel at • 

• E 

i 2p 

AT 



71=0 



n=0 

71=0 




N 



N 

N 




k 

E 

- 71=0 

e rl 

AT 

E Itel AT ■ 

71=0 

• E 

71=0 

G n 

*^rp 

AT 

> 


and (5.13) is used to calculate the linearization error bound for a MIMO system. 


(5.13) 


The error analysis described above is based on the error in the unit impulse response. However, in general, 
the input for the nonlinear model is a step not an impulse; therefore, the impulse error must be derived from 
step response data. In the block diagram shown below, the setup for the error analysis is illustrated. 


5 (t) 



This block diagram can be rearranged so that the outputs of the transfer functions are the impulse responses. 
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Now, the step error is represented as the integration of the impulse error. Therefore, the impulse response 
error can be determined from the step response data by taking the derivative of the step response error. 

e M = J t C V Wirue ~ V (*)) (5-14) 

If the time derivative in the previous expression is approximated by a finite difference expression, then 
equation (5.14) can be rewritten as: 

„ ( (Vtrue ~ Vtrul) ~ fo" ~ V"^) 

B { AT 

and the linearization error bound can be calculated using (5.13) . 

5.2 Model Reduction Error 

When using the modified square root model reduction described in the previous section, the bound on 
the error between the linear full order model and the linear reduced order model can be calculated using the 
method of the previous section. 


(5.15) 


|G(jcu) - G k (jcu)|| oo sC Jjg (t) — g k (t)!^ 


(5.16) 


G (ju) represents the transfer function from the split flux linearization process, and G k {jco) represents the 
transfer function from the model reduction process. Following the procedure outlined above, the error bound 
is calculated from: 


IEOH 


loo ~ °"rnax 


r n 


N 


£ left I AT £K 2 |AT 

n~0 

N 


n=0 

EKilAr Ele&IAr 


n=0 


N 


71=0 


N 


E|e?i|Ar EfelAr 


k L n— 0 


n — 0 


E Kp| A T 

71=0 

e tei AT 

r<— 0 


E | e rp| AT 

71=0 


(5.17) 


and e n is defined as: 


{y n -y n x ) - {y£-y£ l ) 

AT 


(5.18) 
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5.3 Total Error Bound 


As the infinity norm bound of a system is a Nyquist plane bound, the bounds serve as magnitudes of 
error incurred in each step in the modeling process. Thus, the total infinity norm bound is determined as 
the siun of the individual bounds from linearization and model reduction. This is represented as: 

1 1 Gtrue (jv) — Gk O’ w )lloo ^ Untrue (. 7 W ) — G (j w )lloo + !|G(?D) “ Gk (j' w )lloo (5.19) 

where Gtme (jw) represents the transfer function of an exact linear model, G (jco) is the transfer function 
from the linearization process, and Gk (ju>) is the transfer function from the model reduction process. 

6 Example Application 

In this section, two dimensional and three dimensional models of the variable diameter centerbody (VDC) 
inlet are used to illustrate the CFD-based linear modeling technique discussed in this paper. The VDC inlet 
is an axisymmetric inlet with a collapsible centerbody [22] and is typical of inlets being investigated for 
use on aircraft being developed under NASA’s High Speed Civil Transport program. The inlet is designed 
to operate at the following conditions: altitude — 65,000 ft, p^ = 117.8 -jjj, — 395 °R, M ^ — 2.5, 
^compressor face — -33, and 7 — 1.4, and simulation data for the nonlinear model was obtained from the 
CFD code PARC [7]. The steady-state PARC data before the perturbation occurs is used to generate the 
linear model, and the linear model is validated from the transient PARC data that is generated by the 
perturbation. In general, the perturbation is a step input that is not too large so that the small perturbation 
analysis remains valid. Because the CFD model has so many equations (7, 896 for the 2D model and 61, 800 
for the 3D model) it is not numerically feasible to use the entire grid to derive a linear model at this time. 
Therefore, the linear models are calculated based on reduced grids obtained from the original CFD grid. 
These reduced grids are obtained by selecting grid points that approximate the entire grid. Once the linear 
models are generated, the reduced order models are obtained using the modified square root method. The 
nondimensional results for these examples are discussed in the next section, and the data for these matrices 
are listed in section 8.3 of the appendix. 

6.1 Results 

6.1.1 2D VDC Inlet Model: Downstream Mach Number Perturbation 

The two dimensional nonlinear CFD model from cowl lip to compressor face is represented by a 94 x 21 
grid which has 1974 grid points (7896 states). The system has one input and two outputs. The inlet is 
perturbed at 0.01 seconds with a —3% step in compressor face Mach number, and the system outputs for 
static pressure are located downstream of the normal shock at X/R c = 4.08 and X/R c — 5.01. The —3% 
change in the compressor face Mach number produces a 4.24% change in steady-state for the static pressure 
output at XfR c ~ 4.08 and a 3.05% change in the steady-state for the static pressure output at X/R c = 5.01. 
The overall error bounds for various linear models based on reduced grids of the CFD model are listed below, 


Grid 

# FOM States 

# ROM States 

Error Bound 





26 x 3 

312 

13 

0.43461 

26 x 5 

520 

13 

0.43382 

26x7 

728 

13 

0.43342 

48 x 3 

576 

13 

0.41472 

48x5 

960 

13 

0.41449 

48x7 

1344 

13 

0.41424 
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and the nondimensional results for the above models are shown in Figures 1 — 12. 

In Figure 1, the PARC transient data from the 94 x 21 grid is compared to a ROM of the linearized 
CFD model developed from a 26 x 3 grid. In Figure 2, the locations of the FOM eigenvalues are compared 
with those of the ROM. The FOM has 312 states and the ROM has 13 states; this represents a reduction 
of 95.8% in the size of the linear model. The minimum eigenvalue for the FOM is located at -8.9721 x 10 1 
and at —8.9724 x 10 1 for the ROM. The maximum eigenvalue for the FOM is located at —4.6705 x 10 4 and 
at —6.7659 x 10 3 ± j'7.3296 x 10 3 for the ROM. Note the minimum eigenvalue for the FOM and ROM are 
approximately equal, but the maximum eigenvalue for the FOM is further to the left of the imaginary axis 
than for the ROM. This is evident in all of the examples. 

In Figure 3, the PARC transient data from the 94 x 21 grid is compared to a ROM of the linearized 
CFD model developed from a 26 x 5 grid. In Figure 4, the locations of the FOM eigenvalues are compared 
with those of the ROM. The FOM has 520 states and the ROM has 13 states; this represents a reduction 
of 97.5% in the size of the linear model. The minimum eigenvalue for the FOM is located at —8.8813 x 10 1 
and at —8.8783 x 10 1 for the ROM. The maximum eigenvalue for the FOM is located at —87577 x 10 4 and 
at -3.1972 x 10 4 for the ROM. 

In Figure 5, the PARC transient data from the 94 x 21 grid is compared to a ROM of the linearized 
CFD model developed from a 26 x 7 grid. In Figure 6, the locations of the FOM eigenvalues are compared 
with those of the ROM. The FOM has 728 states and the ROM has 13 states; this represents a reduction 
of 98.2% in the size of the linear model. The minimum eigenvalue for the FOM is located at —8.8463 x 10 1 
and at —8.8442 x 10 1 for the ROM. The maximum eigenvalue for the FOM is located at —1.2967 x 10 5 and 
at -2.9908 x 10 4 for the ROM. 

In Figure 7, the PARC transient data from the 94 x 21 grid is compared to a ROM of the linearized 
CFD model developed from a 48 x 3 grid. In Figure 8, the locations of the FOM eigenvalues are compared 
with those of the ROM. The FOM has 576 states and the ROM has 13 states; this represents a reduction 
of 97.7% in the size of the linear model. The minimum eigenvalue for the FOM is located at —1.0385 x 10 2 
and at —1.0398 x 10 2 for the ROM. The maximum eigenvalue for the FOM is located at —8.6676 x 10 4 and 
at -3.9629 x 10 3 ± j8.4573 x 10 3 for the ROM. 

In Figure 9, the PARC transient data from the 94 x 21 grid is compared to a ROM of the linearized 
CFD model developed from a 48 x 5 grid. In Figure 10, the locations of the FOM eigenvalues are compared 
with those of the ROM. The FOM has 960 states and the ROM has 13 states; this represents a reduction 
of 98.6% in the size of the linear model. The minimum eigenvalue for the FOM is located at —1.0399 x 10 2 
and at —1.0424 x 10 2 for the ROM. The maximum eigenvalue for the FOM is located at —1.1310 x 10 5 and 
at -3.4176 x 10 3 ± j8.3316 x 10 3 for the ROM. 

In Figure 11, the PARC transient data from the 94 x 21 grid is compared to a ROM of the linearized 
CFD model developed from a 48 x 7 grid. In Figure 12, the locations of the FOM eigenvalues are compared 
with those of the ROM. The FOM has 1344 states and the ROM has 13 states; this represents a reduction 
of 99% in the size of the linear model. The minimum eigenvalue for the FOM is located at —1.0414 x 10 2 
and at —1.0463 x 10 2 for the ROM. The maximum eigenvalue for the FOM is located at —1.5374 x 10 5 and 
at -3.2790 x 10 3 ± j8.2064 x 10 3 for the ROM. 

6.1.2 3D VDC Inlet Model: Downstream Mach Number Perturbation 

The three dimensional nonlinear CFD model from cowl lip to compressor face is represented by a 103 x 
15 x 8 grid which has 12,360 grid points (61,800 states). The system has one input and two outputs. The 
inlet is perturbed at 0.01 seconds with a —3% step in compressor face Mach number, and the system outputs 
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for static pressure are located at X/R c = 3.9448 and X/R c = 4.9262. The -3% change in the compressor 
face Mach number produces a 3.988% change in steady-state for the static pressure output at X/R c — 3.9448 
and a 3.3362% change in the steady-state for the static pressure output at X/R c = 4.9262. The overall error 
bounds are listed below for two linear models based on reduced grids of the CFD model, 


Grid 

# FOM States 

# ROM States 

Error Bound 





15 x 3 x 3 

675 

17 

0.69196 

20 x 3 x 3 

900 

17 

0.67617 


and the nondimensional results for the above models are shown in Figures 13 — 16. 

In Figure 13, the PARC transient data from the 103 x 15 x 8 grid is compared to a ROM of the linearized 
CFD model developed from a 15 x 3 x 3 grid. In Figure 14, the locations of the FOM eigenvalues are 
compared with those of the ROM. The FOM has 675 states and the ROM has 17 states; this represents 
a reduction of 97.5% in the size of the linear model. The minimum eigenvalue for the FOM is located at 
—1.4375 x 10 2 and at —1.1705 x 10 2 for the ROM. The maximum eigenvalue for the FOM is located at 
-5.4626 x 10 4 and at -8.3988 x 10 3 ± j3.6576 x 10 3 for the ROM. 

In Figure 15, the PARC transient data from the 103 x 15 x 8 grid is compared to a ROM of the linearized 
CFD model developed from a 20 x 3x3 grid. In Figure 16, the locations of the FOM eigenvalues are 
compared with those of the ROM. The FOM has 900 states and the ROM has 17 states; this represents 
a reduction of 98.1% in the size of the linear model. The minimum eigenvalue for the FOM is located at 
—1.1268 x 10 2 and at —1.1268 x 10 2 for the ROM. The maximum eigenvalue for the FOM is located at 
-5.4626 x 10 4 and at -8.4927 x 10 3 ± j3.3323 x 10 3 for the ROM. 

6;2 Discussion 

The CFD steady-state data from the reduced CFD grids seem to provide enough information to calculate 
reasonable linear models that represent the original CFD model. Note that in the two dimensional results, 
an increase in the number of nodes in the y-direction did not have as much effect on lowering the linearization 
error bound as did an increase in the number of nodes in the x-direction. This result is shown in Figure 17 
and illustrates that the inlet model is not strongly two dimensional. In general as the number of grid points is 
increased in the reduced CFD grid, the better the linear model; in other words, linear models obtained from 
larger reduced CFD grids have lower linearization error bounds. For the three dimensional model, increasing 
the number of nodes in the x-direction lowered the linearization error bound but not by very much. This is 
probably due to the fact that so little of the original CFD grid is being used in the development of the linear 
model. The three dimensional example is included mainly as an application of the linear modeling concept 
and model reduction process. 

Linear results based on one dimensional data obtained from averaged two dimensional data are shown 
in Figures 18 and 19. As was noted above, the inlet model is not strongly two dimensional; therefore, the 
dynamic behavior of the inlet can be represented by a one dimensional linear model. The linear model 
based on the averaged steady-state data has 282 states. The linear results based on one dimensional data 
obtained from averaged three dimensional data are shown in Figures 20 and 21. The linear model based 
on the averaged steady-state data has 309 states. Again, the three dimensional inlet model is adequately 
represented by a one dimensional linear model based on the averaged steady-state data. 

As was noted in the one dimensional modeling paper [9], the full order model eigenvalues closest to 
the origin are reproduced in the reduced order model, but the eigenvalues farther away from the origin 
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are approximated by just a few eigenvalues in the reduced order model. The model reduction for the 
linearized models based on the multidimensional CFD data is not as straight forward as it is for the one 
dimensional models. In the one dimensional case, the model reduction process for the linear models used 
established methods. For the multidimensional linear models, the model reduction process had to be modified 
by implementing approximation methods for the Grammian and SVD calculations. These approximation 
methods require that the user declare error tolerances for the algorithms and approximation sizes for the 
matrices. The smaller the error tolerance used, the longer the approximation method takes to converge; so 
there is a trade-off between computation speed and accuracy. In addition, as the full order linear models 
become larger, the approximation sizes need to be made larger; this also increases the amount of time needed 
to compute the reduced order model. Unfortunately, there are no proven guidelines that can be used when 
choosing the approximation sizes and error tolerances for the model reduction method; improvement in this 
process needs to be addressed. 

7 Conclusion 

Linear dynamic models of multidimensional mixed compression inlets have been developed from steady- 
state CFD results. As was the case for the one dimensional CFD models, it is possible to obtain a linear model 
from the spatial information and steady-state operating conditions of a multidimensional CFD simulation. 
The linearization process for the multidimensional models was more difficult due to the different types of 
boundary conditions. The small perturbation models that result are useful for control applications, real-time 
simulations, and model reduction. Additional types of inputs need to be developed as was done for the linear 
models based on quasi-one dimensional CFD results [9]. 

Model reduction proved to be a more difficult problem than for the one dimensional case due to the 
size of the computational grids of the multidimensional CFD models. Modifications to the existing model 
reduction method were made which produced viable reduced order models. Furthermore, model reduction 
for systems on the order of a thousand states is possible. Work remains to be completed on more efficient 
mbdel reduction algorithms for these large scale systems. 
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8 Appendices 


8.1 Symbols 


A 

A 

B 

B 

Bbc 

C 

C 

D 

c 

E 

E 

e 

e(t) 

e 

7 

G 

~9 

sW 

h 

K 

3 

M 

m 

N 2 

n 3 

V 

p 

Q 

R c 

s 

t 


w 

~u 

w 

X 

y 

y(t) 


Jacobian for 
system matrix 
Jacobian for 
input matrix 

physical boundary condition 

Jacobian for 
output matrix 
input/output matrix 

speed of sound 

total energy per unit mass 
error transfer function matrix 
internal energy per unit mass 


error response 

error impulse matrix 

flux vector for the x-direction 

transfer function 

flux vector for the y-direction 

impulse response 

flux vector for the z-direction 

matrix of right eigenvectors 

\/-T 

Mach number 

massflow per unit area (^rfi 2 j 
number of grid points in the x-direction 
number of grid points in the y-direction 
number of grid points in the z-direction 

pressure 

controllability Grammian 
observability Grammian 
cowl radius (ft) 
entropy 
time (sec) 

velocity ( ^ j in the x-direction 
velocity ( ^ ) in the y-direction 

velocity ( ) in the z-direction 

vector of conservative variables 
velocity vector 

vector of characteristic variables 

distance (ft) 

distance (ft) 

output reponse 

distance (ft) 
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8.1.1 Greek Symbols 


AT sampling time (sec) 
e e — pE, total energy per unit volume 

8 small perturbation variable 
8 ( t ) unit impulse 

7 ratio of specific heats 

A diagonal matrix of local characteristics 

A eigenvalue 

p density 

a singular value 

lu ■ frequency (^ 7 ) 

8.1.2 Subscripts 

be boundary condition 
i grid point in the x- direction 

j grid point in the y-direction 
k grid point in the z-direction 
k size of reduced order model 
L left 

p number of outputs 
R right 

r number of inputs 
t total conditions 

8.1.3 Superscripts 

N numerical 
n time level 
P physical 
T transpose 



8.2 Figures 


8.2.1 2D VDC Inlet Model: Downstream Mach Number Perturbation 
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8.2.2 3D VDC Inlet Model: Downstream Mach Number Perturbation 
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8.2.3 2D VDC Inlet Model: Data Comparison 
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600 800 

Number of States 


Figure 17: 2D VDC Inlet Linear Model Comparison 
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8.2.4 2D VDC Inlet Model: ID Model from 2D Averaged Data 
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8.2.5 3D VDC Inlet Model: ID Model from 3D Averaged Data 
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8.3 Matrices for Reduced Order Linear Models 

8.3.1 2D VDC Inlet Model: Downstream Mach Number Perturbation 

8.3.1. 1 26 x 3 Grid 
System Matrix [A] 


Columns 1 through 6 
1 . 0255e+002 


2 . 2093e+002 
•1 . 8683e+002 
1 . 8411e+002 
1 . 3245e+002 
1 . 0913e+002 
2 . 2540e+001 
4 . 043 8e+001 
2 . 17 64e+001 
2 . 2328e+001 
1 . 5298e+001 
4 . 0977e+000 
2 . 5383e+000 


-1 . 4033e+002 
-3 . 7917e+002 
1 . 4195e+003 
-9 . 2390e+002 
4 . 8332e+002 
-4 . 4412e+002 
-9 . 1965e+001 
1. 6202e+002 
-8 . 7417e+001 
-8 . 973 le+001 
6 . 1472e+001 
-1 . 6449e+001 
1 . 0197e+001 


1 . 3501e+002 
-1. 1727e+003 
-4 . 6557e+002 
7 . 1603e+002 
-8 . 143 8e+002 
6 . 3725e+002 
1 . 2418e+002 
-2 . 2395e+002 
1 . 2041e+002 
1 . 2412e+002 
-8 . 5217e+001 
2 . 2830e+001 
-1 . 4137e+001 


Columns 7 through 12 


1. 0069e+001 
1 . 1725e+002 
8 . 9505e+001 
3 . 4241e+002 
1 . 0437e+003 
7 . 7816e+002 
4 . 8341e+002 
2 . 5063e+003 
1. 0047e+003 
9 . 03 63e+002 
6 . 1240e+002 
1 . 7540e+002 
1 . 0903e+002 


4 . 0819e+001 
5 . 8929e+001 
-1 . 7056e+002 
3 . 0550e+002 
-6 . 8509e+001 
-3 . 7345e+003 
-1. 9894e+003 
-3 . 9720e+003 
4 . 1477e+003 
3 . 6098e+003 
-3 . 1572e+003 
8 . 1722e+002 
-5 . 042 6e+002 


-1 . 9063e+001 
1. 6378e+001 
1 . 4244e+002 
-3 . 4980e+002 
-4 . 0756e+002 
2 . 0909e+003 
2 . 3277e+002 
1. 5434e+002 
-1 . 413 9e+003 
-1. 6716e+003 
1 . 8972e+003 
-5 . 2263e+002 
3 . 4477e+002 


Column 13 
4 . 4124e-001 
7 . 9507e+000 
9 . 3206e+000 
-2 . 8351e+001 
-5 . 6142e+001 
-8 . 3960e+000 
-2 . 5088e+002 
1. 6218e+002 
-3 . 4820e+002 
-1. 9747e+003 
1 . 7715e+003 
-1. 9314e+003 
-5 . 4131e+003 


-8 . 3608e+001 
7 . 5097e+002 
6 . 5513e+002 
-1 . 2450e+003 
2 . 6420e+003 
-1. 6596e+003 
-3 . 0009e+002 
5 . 5744e+002 
-2 . 9822e+002 
-3 . 0930e+002 
2 . 1272e+002 
-5 . 7033e+001 
3 . 53 68e+001 


-1 . 9510e+001 
6 . 1803e+001 
1. 6331e+002 
-3 . 7484e+002 
-5 . 4782e+002 
2 . 2577e+003 
-5. 6317e+002 
3 . 2215e+003 
-4 . 1516e+003 
-5. 0254e+003 
8 . 1878e+003 
-1 . 9743e+003 
1 . 2017e+003 


3 . 4645e+001 
4 . 9048e+002 
5 . 5007e+002 
-2 . 3700e+003 
-1 . 9782e+003 
3 . 5833e+003 
8 . 0327e+002 
-1 . 2103e+003 
6 . 6887e+002 
6 . 7680e+002 
-4 . 6229e+002 
1 . 2324e+002 
-7. 6172e+001 


-7 . 483 le+000 
-6 . 7737e+001 
-2 . 7573e+001 
1. 1328e+002 
4 . 3561e+002 
2 . 3524e+002 
8 . 7657e+002 
1 . 2781e+003 
4 . 1537e+002 
-2 . 6730e+003 
-4 . 3154e+003 
2 . 0910e+003 
-1. 6925e+003 


1 . 0379e+002 
1 . 0045e+002 
-5 . 5065e+002 
1 . 2228e+003 
-9 . 0652e+002 
-4 . 3863e+003 
-2 . 1006e+003 
3 . 7847e+003 
-1. 9898e+003 
-1. 8642e+003 
1.2176e+003 
-3 . 2653e+002 
2 . 0254e+002 


3 . 6329e+000 
5 . 8024e+000 
- 1 . 3 655e+001 
1 . 8785e+001 
-3 . 0598e+001 
-2 . 4656e+002 
1 . 5314e+001 
-9 . 7361e+002 
5 . 7451e+002 
2 . 2351e+003 
1 . 9289e+002 
-1 . 5182e+003 
3 . 5735e+003 
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Input Matrix [B] 

9 . 6540e+002 
-8.8858e+002 
7 . 8651e+002 
-8 . 419 6e+002 
6 . 1995e+002 
-5 . 1694e+002 
-1.0622e+002 
1 . 9015e+002 
-1. 0237e+002 
-1.0507e+002 
7 . 2004e+001 
-1. 9283e+001 
1 . 1952e+001 


Output Matrix [C] 

Columns 1 through 6 

-6.77 67e+002 -8.6623e+002 -7.4260e+001 3.2476e+002 5.3850e+002 

-6 . 8859e+002 2.7100e+002 7.8304e+002 -7.8850e+002 -3.1708e+002 

Columns 7 through 12 

1 . 2966e+002 1.1676e+002 -7.9688e+000 1.3084e+001 -6.9400e+001 

-6 . 2663e+001 1.5445e+002 -1.2043e+002 -1.3748e+002 1.6733e+001 

Column 13 
9 . 9737e+000 
-6 . 7083e+000 


Input/Output Matrix [D] 

o 

0 


2 . 3275e+002 
4 . 6508e+002 


1 . 1069e+001 
1 . 2464e+001 
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8.3.1. 2 26 x 5 Grid 


System Matrix [A] 


Columns 1 through 6 


-1 . 013 0e+002 
2 . 162 0e+002 
-1 . 8431e+002 
1. 8193e+002 
1 . 307 5e+002 
-1 . 0768e+002 
2 . 6730e+001 
-4 . 4093e+001 
9 . 7886e+000 
-2 . 1988e+001 
-1 . 3105e+001 
-1 . 7192e+000 
4 . 1066e+000 


-1 . 3867e+002 
-3 . 7095e+002 
1 . 4165e+003 
-9 . 1931e+002 
-4 . 7567e+002 
4 . 3814e+002 
-1 . 0868e+002 
1 . 7667e+002 
-3 . 9404e+001 
8 . 8355e+001 
5 . 263 6e+001 
6 . 9168e+000 
-1. 6513e+001 


1 . 3372e+002 
-1 . 1746e+003 
-4 . 5959e+002 
7 . 0664e+002 
8 . 0371e+002 
-6 . 2956e+002 
1 . 4745e+002 
-2 . 4447e+002 
5 . 4141e+001 
-1 . 2237e+002 
-7 . 3084e+001 
-9 . 5971e+000 
2 . 2911e+001 


-8 . 4302e+001 
7. 6559e+002 
6 . 5398e+002 
-1 . 2391e+003 
-2. 6390e+003 
1 . 6548e+003 
-3 . 6009e+002 
6 . 1256e+002 
-1 . 3411e+002 
3 . 0707e+002 
1 . 8373e+002 
2 . 4089e+001 
-5 . 7641e+001 


-3 . 5668e+001 - 
-4.9205e+002 - 
-5 . 3 612e+002 
2 . 3 686e+003 - 
-1 . 9474e+003 - 
3 . 5253e+003 - 
-9 . 2160e+002 
1 . 3 177e+003 . - 
-3 . 0788e+002 
6 . 6459e+002 - 
3 . 9430e+002 - 
5 . 1498e+001 - 
-1 . 2292e+002 


Columns 7 through 12 


7 . 2251e+000 
1 . 1224e+002 
1 . 0041e+002 
3 . 7172e+002 
1 . 0192e+003 
1 . 0833e+003 
6 . 9951e+002 
2 . 8122e+003 
4 . 1704e+002 
1. 0366e+003 
6 . 5647e+002 
9 . 0397e+001 
2 . 163 3e+002 


4 . 3 628e+001 
6 . 0941e+001 
-1 . 8714e+002 
3 . 4665e+002 
5 . 0733e+001 
4 . 0840e+003 
-1 . 349 6e+003 
-5 . 0362e+003 
3 . 3567e+003 
-4 . 0932e+003 
-3 . 1621e+003 
-3 . 9964e+002 
9. 5083e+002 


-4 . 6580e+000 
3 . 2746e+001 
7 . 9812e+001 
-2 . 4685e+002 
3 . 9817e+002 
-7 . 5038e+002 
-3 . 0139e+002 
-1 . 3600e+003 
-2 . 92 82e+002 
3 . 3911e+002 
6 . 7849e+002 
1 . 1878e+002 
-2 . 8247e+002 


1. 8437e+001 
-7 . 8357e+001 
-1. 6431e+002 
3 . 7983e+002 
-5 . 9550e+002 
2 . 2990e+003 
1 . 3548e+003 
-4. 1973e+003 
3 . 7232e+003 
-5. 6575e+003 
-8 . 0595e+003 
-9 . 1242e+002 
2 . 2335e+003 


-9 . 0247e+000 - 
-5 . 4084e+001 - 
-4 . 9668e+000 
5 . 0626e+001 - 
-3 . 2707e+002 - 
-4 . 5683e+002 - 
7 . 3993e+002 - 
1 . 8419e+003 
7 . 4581e+002 - 
3 . 2707e+003 
-4 . 2 062e+003 - 
-1 . 3702e+003 - 
3 . 5142e+003 


Column 13 
3 . 4407e+000 
-1.967 6e+001 
-5 . 0931e+001 
1 . 0651e+002 
-1 . 5241e+002 
2 . 8627e+002 
5 . 9106e+002 
-1 . 7558e+003 
9 . 0314e+002 
-9 . 9978e+003 
-3 . 3404e+003 
7 . 7656e+003 
-3 . 2657e+004 


1 . 0265e+002 
9 . 9622e+001 
5 . 522 le+002 
1 . 2605e+003 
8 . 9620e+002 
4 . 3559e+003 
2 . 5613e+003 
4 . 1251e+003 
8 . 7511e+002 
1 . 8396e+003 
1 . 0467e+003 
1 . 3869e+002 
3 . 3069e+002 


2 , 2976e+000 
7 . 3942e-001 
1 . 0593e+001 
8 . 9564e+000 
1 . 003 9e+001 
1. 6212e+002 
8 . 9542e+001 
7. 1874e+002 
3 . 0943e+002 
1 . 2195e+003 
8. 6596e+002 
1 . 6653e+003 
4 . 6800e+003 


58 



Input Matrix [B] 

9 . 5925e+002 
-8.7457e+002 
7 . 8027e+002 
-8 . 3652e+002 
-6 . 1544e+002 
5 . 1305e+002 
-1 . 2664e+002 
2 . 0856e+002 
-4 . 6323e+001 
1 . 0408e+002 
6 . 2041e+001 
8 . 1399e+000 
-1. 9454e+001 


Output Matrix [C] 

Columns 1 through 6 

-6 . 7501e+002 -8.6288e+002 -6.8448e+001 3.2103e+002 -5.3706e+002 -2.3012e+002 

-6 . 8408e+002 2.7743e+002 7.7683e+002 -7.8971e+002 3.0941e+002 -4.6601e+002 

Columns 7 through 12 

1 . 2159e+002 1.2135e+002 3.0474e+001 -2.4111e+001 -6.1034e+001 -5.1730e+000 

-7 . 4744e+001 1.7143e+002 -6.4954e+001 1.3897e+002 -1.3402e+000 -9.1019e+000 

Column 13 
-1 . 7028e+001 
3 . 8283e+001 


Input/Output Matrix [D] 

o 

0 
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8.3. 1.3 26 x 7 Grid 


System Matrix [A] 


Columns 1 through 6 
1 . 0094e+002 -1.3798e+002 


2 . 1573e+002 
•1 . 8327e+002 
1 . 8062e+002 
1 . 3007e+002 
1 . 0753e+002 
2 . 7303e+001 
4 . 4105e+001 
5 . 5409e+000 
2 . 3940e+001 
1 . 2842e+001 
1 . 3565e+000 
3 . 9972e+000 


-3 . 7070e+002 
1 . 4143e+003 
-9. 1530e+002 
-4 . 7378e+002 
4 . 3790e+002 
-1 . 1115e+002 
1 . 7 693e+002 
-2 . 2388e+001 
9 . 6304e+001 
5 . 1639e+001 
5 . 4588e+000 
-1 . 6089e+001 


1 . 3332e+002 
-1 . 1756e+003 
-4 . 5510e+002 
6 . 9857e+002 
7 . 9585e+002 
-6.2659e+002 
1 . 4999e+002 
-2 . 4365e+002 
3 . 0435e+001 
-1 . 3274e+002 
-7 . 1353e+001 
-7 . 5427e+000 
2 . 2219e+001 


-8 . 4070e+001 
7. 6819e+002 
6 . 4886e+002 
-1 . 2265e+003 
-2 . 6273e+003 
1. 6508e+003 
-3 . 6590e+002 
6. 1090e+002 
-7 . 4701e+001 
3 . 3332e+002 
1 . 7951e+002 
1 . 8946e+001 
-5 . 5950e+001 


-3 . 6817e+001 
-4 . 8 668e+002 
-5 . 2640e+002 
2 . 3 685e+003 
-1 . 9263e+003 
3 . 5024e+003 
-9 . 3486e+002 
1 . 3 110e+003 
-1 . 7837e+002 
7 . 1945e+002 
3 . 8441e+002 
4 . 0374e+001 
-1 . 1891e+002 


Columns 7 t. 
7 . 3147e+000 
1 . 1409e+002 
1 . 0414e+002 
-3 . 9134e+002 
1 . 0075e+003 
-1 . 1153e+003 
-7 . 1866e+002 
2 . 7725e+003 
-1. 8722e+002 
1 . 1408e+003 
6 . 4883e+002 
7 . 1734e+001 
-2 . 117 0e+002 


trough 12 
4 . 3408e+001 
5 . 8557e+001 
-1 . 8876e+002 
3 . 5501e+002 
5 . 1882e+001 
4. 0569e+003 
-1 . 2480e+003 
-5. 1324e+003 
2 . 8397e+003 
-4 . 5755e+003 
-3 . 1493e+003 
-3 . 2147e+002 
9 . 4532e+002 


-1 . 3041e+000 
3 . 3388e+001 
6 . 2528e+001 
-2 . 1222e+002 
3 . 5019e+002 
-3 . 9396e+002 
-3 . 2944e+002 
-1 . 8005e+003 
-8. 9068e+001 
-2 . 1092e+001 
3 . 1203e+002 
4 . 9978e+001 
-1 . 5810e+002 


1 . 7895e+001 
-7 . 1789e+001 
-1 . 6068e+002 
3 . 7902e+002 
-5. 5588e+002 
2 . 2047e+003 
1 . 3015e+003 
-3 . 9944e+003 
2 . 8489e+003 
-6 . 0759e+003 
-7 . 6279e+003 
-6 . 9491e+002 
2 . 1186e+003 


-8 . 8425e+000 - 
-5 . 3 681e+001 - 
-2 . 8 642e+000 
4 . 3260e+001 - 
-3 . 2941e+002 - 
-4 . 2147e+002 - 
7 . 1261e+002 - 
1 . 7808e+003 
8 . 6305e+002 - 
3 . 0742e+003 
-4 . 1982e+003 - 
-9 . 2567e+002 - 
3 . 3551e+003 


Column 13 
3 . 7291e+000 
-1 . 9685e+001 
-4 . 8379e+001 
9 . 923 8e+001 
-2 . 0361e+002 
4 . 3820e+002 
5 . 2903e+002 
-1 . 9198e+003 
6 . 2660e+002 
-8 . 8375e+003 
-3 . 7006e+003 
1. 9702e+003 
-3 . 1368e+004 


1 . 0216e+002 
9 . 5743e+001 
5 . 5463e+002 
1 . 2827e+003 
8.9618e+002 
4 . 3563e+003 
2 . 6306e+003 
4. 1186e+003 
4. 9365e+002 
2 . 0037e+003 
1 . 0273e+003 
1 . 095 6e+002 
3 . 22 62e+002 


2 . 4964e+000 
1 . 9640e+000 
9 . 1749e+000 
2 . 8165e+000 
5 . 1386e+000 
2 . 0042e+002 
6 . 5397e+001 
7 . 9735e+002 
2 . 2395e+002 
1. 1337e+003 
9 . 4632e+002 
1. 0687e+003 
5 . 1491e+003 
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Input Matrix [B] 

9 . 573 8e+002 
-8 . 7445e+002 
7 . 7681e+002 
-8 . 3179e+002 
-6 . 1311e+002 
5 . 1318e+002 
-1 . 2957e+002 
2 . 0897e+002 
-2 . 6269e+001 
1. 1351e+002 
6 . 0896e+001 
6 . 4309e+000 
-1. 8967e+001 


Output Matrix [C] 

Columns 1 through 6 

-6.7403e+002 -8.6047e+002 -6.7453e+001 3.2312e+002 -5.3624e+002 -2.2718e+002 

-6 . 8243e+002 2.7756e+002 7.7503e+002 -7.9089e+002 3.0228e+002 -4.6766e+002 

Columns 7 through 12 

1 . 2373e+002 1.1940e+002 3.6240e+001 -2.3144e+001 -5.9665e+001 -6.7813e+000 

-7 . 7027e+001 1.7281e+002 -4.8324e+001 1.3570e+002 -1.8189e+000 -9.3922e+000 

Column 13 
-1 . 5488e+001 
4 . 1400e+001 


Input/Output Matrix [D] 

0 

0 
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8.3.1.4 48 x 3 Grid 


System Matrix [A] 


Columns 1 through 6 

■1 . 0483e+002 -1.9628e+002 


2 . 5622e+002 
2 . 3583e+002 
1.5793e+002 
1 . 6230e+002 
1 . 3473e+002 
5 . 8742e+000 
8 . 1706e+001 
5 . 0294e+001 
1 . 8189e+001 
2 . 2756e+001 
1 . 4650e+001 
6 . 4721e+000 


-2 . 8724e+002 
1 . 4265e+003 
-5 . 2051e+002 
4 . 0304e+002 
-4 . 2055e+002 
1. 8187e+001 
2 . 3803e+002 
1 . 4895e+002 
5 . 3657e+001 
6 . 7384e+001 
-4 . 3277e+001 
-1 . 9292e+001 


1 . 2608e+002 
-1 . 2499e+003 
-5 . 2074e+002 
5 . 1647e+002 
-1 . 0421e+003 
7 . 1668e+002 
-2 . 4751e+001 
-3 . 9587e+002 
-2 . 4552e+002 
-9 . 0158e+001 
-1 . 1268e+002 
7 . 2654e+001 
3 . 2306e+001 


-1. 0084e+002 
2 . 7453e+002 
4 . 7068e+002 
-5 . 5848e+002 
2 . 3818e+003 
-1 . 2319e+003 
3 . 2393e+001 
5 . 8413e+002 
3 . 6853e+002 
1 . 3775e+002 
1 . 7121e+002 
-1 . 1063e+002 
-4. 8944e+001 


3 . 2454e+001 
4 . 4786e+002 
9 . 9519e+002 - 
-2 . 2011e+003 
-1 . 2281e+003 - 
2 . 5626e+003 - 
-1 . 4856e+002 
-1 . 2593e+003 
-8 . 0614e+002 
-2 . 8291e+002 
-3 . 4955e+002 
2 . 2303e+002 - 
9 . 8268e+001 - 


Columns 7 through 12 


9 . 6708e+000 
-5 . 5568e+001 
-1 . 4053e+002 
1 . 8357e+002 
5 . 1149e+002 
-7 . 9759e+002 
-1 . 4526e+001 
-2 . 1361e+003 
3 . 827 6e+001 
-1 . 392 8e+002 
-1 . 2722e+002 
8 . 1512e+001 
1 . 1508e+001 


8 . 8381e+001 
1 . 0082e+002 
-2 . 2539e+002 
4 . 4622e+002 
1. 5517e+002 
-4 . 2275e+003 
1 . 823 6e+003 
-2 . 7675e+003 
-3 . 1697e+003 
-1 . 5628e+003 
-1. 9778e+003 
1 . 0234e+003 
4 . 4335e+002 


2 . 8703e+001 
-1 . 0002e+002 
-2 . 9719e+002 
3 . 9332e+002 
8 . 7505e+002 
-1 . 5345e+003 
-7 . 6998e+002 
-6 . 5346e+002 
-2 . 7704e+003 
-3 . 5625e+003 
-3 . 6405e+003 
1. 8274e+003 
7 . 3170e+002 


-1 . 7 045e+001 
-9 . 7755e+000 
5. 8658e+001 
-1. 2578e+002 
8 . 8372e+001 
4 . 5363e+002 
1. 6516e+002 
1 . 4433e+003 
2 . 8469e+003 
-2 . 8863e+002 
1 . 5949e+003 
9 . 8400e+002 
8 . 6116e+002 


-3 . 6775e+001 
-9 . 4945e+001 
1 . 9830e+001 
-5.687 5e+001 
1 . 5213e+000 
1 . 2998e+003 
4 . 53 97e+001 
3 . 1795e+003 
3 . 3216e+003 
-4 . 5229e+003 
-1 . 4218e+003 
2 . 2969e+003 
1 . 2 048e+003 


Column 13 
-3 . 2506e+000 
1. 0952e+001 
6 . 8427e+001 
-1 . 2273e+002 
-1 . 4504e+002 
-1 . 1634e+001 
6 . 2712e+002 
8 . 9129e+002 
3 . 0751e+003 
-1 . 3937e+003 
1 . 1527e+003 
5 . 3541e+003 
-4 . 2265e+002 


1 . 33 9 9e+002 
2 . 3388e+002 
3 . 6918e+002 
9 . 3637e+002 
1 . 1625e+003 
2 . 2635e+003 
4 . 2708e+002 
3 . 8865e+003 
1 . 8197e+003 
5 . 6571e+002 
6 . 9668e+002 
4 . 8167e+002 
2 , 1911e+002 


1 . 7547e+001 
3 . 7994e+001 
1 . 8108e+002 
2 . 7 848e+002 
5 . 9753e+002 
5 . 3892e+002 
1 . 5606e+003 
2 . 7789e+003 
7 . 1778e+003 
2 . 9258e+003 
1. 5476e+003 
9 . 1383e+003 
8 . 5998e+003 
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Input Matrix [B] 

9 . 9631e+002 
-9 . 2447e+002 
9 . 9931e+002 
-7 . 053 6e+002 
7. 6353e+002 
-6 . 5291e+002 
2 . 7642e+001 
3 . 8592e+002 
2 . 3861e+002 
8 . 6305e+001 
1 . 0810e+002 
-6 . 9656e+001 
-3 . 0899e+001 


Output Matrix [C] 

Columns 1 through 6 

-6 . 8180e+002 -9.3312e+002 -4.1095e+002 1.1514e+002 6.9145e+002 

-7 . 2794e+002 1.3294e+002 8.9597e+002 -6.7167e+002 -4.7058e+002 

Columns 7 through 12 

-7 . 5019e+001 2.2344e+002 -1.0630e+002 -4 . 2305e+001 -1.4597e+002 - 

1 . 3009e+002 3.5647e+002 2.7348e+002 -6.5985e+001 -1.0972e+002 

Column 13 
2 . 6625e+001 
-5 . 1848e+001 


Input/Output Matrix [D] 

0 

0 


3 . 9987e+002 
5 . 2363e+002 


4 . 4552e+001 
1 . 6872e+002 
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8. 3.1. 5 48 x 5 Grid 


System Matrix [A] 


Columns 1 through 6 
1 . 0455e+002 -1.9623e+002 
-2 . 7304e+002 


2 . 4939e+002 
2 . 3701e+002 
1 . 6106e+002 
1 . 6248e+002 
1 . 3445e+002 
2 . 1338e+001 
7 . 5448e+001 
4 . 6193e+001 
3 . 6405e+001 
2 . 1556e+000 
1 . 3818e+001 
7 . 3590e+000 


1 . 4268e+003 
-5 . 1776e+002 
3 . 8723e+002 
-4 . 1173e+002 
5 . 9735e+001 
-2 . 1514e+002 
1 . 3427e+002 
-1 . 0590e+002 
6. 5082e+000 
-4 . 003 8e+001 
2 . 1405e+001 


1 . 2796e+002 
-1 . 2457e+003 
-5 . 2322e+002 
5 . 3116e+002 
-1 . 0488e+003 
7 . 1160e+002 
-1 . 0788e+002 
3 . 6488e+002 
-2 . 2532e+002 
1 . 8048e+002 
-1 . 0495e+001 
6. 8588e+001 
-3 . 6831e+001 


-1 . 0679e+002 
2 . 6547e+002 
4 . 7425e+002 
-5 . 832 8e+002 
2 . 4369e+003 
-1 . 2488e+003 
1. 6919e+002 
-5 . 4693e+002 
3 . 4517e+002 
-2 . 8175e+002 
1. 6038e+001 
-1 . 0712e+002 
5 . 7324e+001 


Columns 7 through 12 

-8 . 8483e+001 
-1 . 2600e+002 


1 . 8045e+001 
9 . 453 0e+001 
8 . 5641e+001 
8 . 1814e+001 
4 . 619 6e+002 
5 . 0774e+002 
1. 8170e+002 
1 . 2323e+003 
8 . 3553e+002 
7 . 5695e+002 
9 . 3351e+001 
2.202 6e+002 
1 . 3483e+002 


2 . 67 65e+002 
-5 . 5609e+002 
-4 . 6918e+002 
4 . 3305e+003 
2 . 7981e+003 
-2 . 4868e+003 
2 . 6807e+003 
-3 . 0490e+003 
9 . 2561e+001 
-9 . 5159e+002 
4 . 5876e+002 


3 . 2974e+001 
-8 . 2547e+001 
-3 . 0785e+002 
4 . 4131e+002 
9 . 0296e+002 
-1 . 7812e+003 
1 . 6803e+002 
9 . 2563e+002 
-2 . 2497e+003 
5 . 2248e+003 
1 . 5412e+002 
1 . 5415e+003 
~7 . 2494e+002 


3 . 4044e+001 
8 . 7920e+001 
-7 . 4962e+000 
8 . 6144e+001 
-2 . 1406e+002 
-1. 1103e+003 
-8. 2331e+002 
2 . 9132e+003 
-3 . 4683e+003 
-2 . 5891e+003 
-2 . 2558e+003 
-2 . 2074e+003 
1. 6582e+003 


Column 13 
5 . 573 0e+000 
-1 . 5162e+001 
-6 . 1750e+001 
9 . 8114e+001 
1 . 0258e+002 
-4 . 2661e+001 
4 . 0652e+002 
1 . 3328e+003 
-3 . 0297e+003 
-5 . 2930e+002 
-1 . 5980e+003 
-4 . 7722e+003 
3 . 4494e+003 


2 . 8222e+001 
4 . 8648e+002 
1 . 0004e+003 - 
-2 . 1959e+003 
-1 . 2197e+003 - 
2 . 5172e+003 - 
-2 . 6683e+002 
1 . 1850e+003 - 
-7 . 3878e+002 
5 . 5866e+002 - 
-3 . 3701e+001 
2 . 0881e+002 - 
-1 . 1005e+002 


-1 . 09 68e+001 
-1 . 1126e+002 - 
-1 . 1002e+002 - 
9 . 9207e+001 
1 . 4286e+002 
6 . 5584e+002 - 
5 . 1974e+002 
2 . 7243e+002 
-1 . 2280e+003 - 
4 . 5829e+003 - 
7 . 2950e+001 - 
4 . 2164e+002 - 
2 . 8144e+002 


1 . 3286e+002 
2 . 2737e+002 
3 . 2782e+002 
9 . 02 65e+002 
1. 1346e+003 
2 . 2774e+003 
8 . 623 5e+002 
3 . 6783e+003 
1 . 7351e+003 
1. 1426e+003 
5 . 1640e+001 
4 . 5851e+002 
2.513 le+002 


1 . 5509e+001 
6 . 8768e+001 
1 . 7561e+002 
2 . 5569e+002 
5 . 1793e+002 
3 . 519 6e+002 
1 . 1448e+003 
3 . 3454e+003 
7 . 4 550e+003 
2 . 7171e+002 
2 . 4987e+003 
8 . 4456e+003 
1 . 0547e+004 
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Input Matrix [B] 

9. 9516e+002 
-8 . 9662e+002 
1 . 0041e+003 
-7 . 1895e+002 
7 . 6645e+002 
-6.5235e+002 
1 . 0104e+002 
-3 . 5682e+002 
2 . 1942e+002 
-1 . 7330e+002 
1 . 0295e+001 
-6 . 5821e+001 
3 . 5145e+001 


Output Matrix [C] 

Columns 1 through 6 

-6 . 8178e+002 -9.4470e+002 -3.9999e+002 9.1059e+001 6.8014e+002 

-7 . 2780e+002 1.4996e+002 8.9168e+002 -6.7764e+002 -4.7990e+002 

Columns 7 through 12 

1. 5453 e+ 002 -2.0383e+002 -8.8999e+001 1.5096e+002 -1.1461e+002 - 

-2 . 7551e+001 -3.8939e+002 2.8982e+002 8.4253e+001 1.0722e+001 

Column 13 
-1. 8861e+001 
5 . 3164e+001 


Input/Output Matrix [D] 

o 

o 


4 . 0995e+002 
5 . 1178e+002 


5 . 8842e+001 
1 . 6015e+002 
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8. 3.1. 6 48 x 7 Grid 


System Matrix [A] 


Columns 1 through 6 

~1.9907e+002 
-2 . 6835e+002 


1. 0409e+002 
2 . 483 6e+002 
2 . 3417e+002 
1 . 5853e+002 
1. 6412e+002 
1 . 3282e+002 
6 . 9811e+001 
5 . 7087e+001 
4 . 7621e+001 
3 . 0830e+001 
1 . 2461e+001 
1. 7659 e+ 001 
3 . 8254e+000 


1 . 4292e+003 
-5 . 1215e+002 
3 . 9156e+002 
-4 . 0580e+002 
1 . 9870e+002 
-1. 6083e+002 
1 . 3787e+002 
8 . 8913e+001 
3 . 6196e+001 
5 . 0840e+001 
-1. 1058e+001 


1 . 2492e+002 
-1 . 2575e+003 
-5 . 0954e+002 
5 . 0959e+002 
-1 . 0466e+003 
6 . 9360e+002 
-3 . 3085e+002 
2 . 7634e+002 
-2 . 2965e+002 
-1 . 5121e+002 
-6 . 0560e+001 
-8. 6416e+001 
1 . 9034e+001 


-1. 1088e+002 
2. 8630e+002 
4 . 7450e+002 
-5. 6875e+002 
2 . 4773e+003 
-1 . 2409e+003 
4. 9540e+002 
-4 . 2284e+002 
3 . 5650e+002 
2 . 3932e+002 
9 . 4811e+001 
1. 3609e+002 
-2 . 9548e+001 


2 . 9894e+001 
4 . 5674e+002 
9 . 9335e+002 - 
-2 . 1924e+003 
-1 . 2193e+003 - 
2 . 4404e+003 - 
~1.1131e+003 
8 . 3062e+002 - 
-7 . 5180e+002 
-4 . 6115e+002 
-1 . 9195e+002 
-2 . 6390e+002 
5 . 6577e+001 - 


Columns 7 through 12 

-5 . 7400e+001 
-1 . 4669e+002 


6 . 5031e+001 
5 . 3521e+001 
2 . 7061e+002 
5 . 4572e+002 
5 . 9517e+002 
3 . 3884e+003 
1. 8312e+003 
3 . 2844e+003 
2 . 0740e+003 
2 . 0697e+003 
5 . 7788e+002 
9 . 4256e+002 
1 . 903 6e+002 


4 . 4545e+001 
-2 . 1265e+002 
1 . 9508e+002 
2 . 4739e+003 
-3 . 2676e+002 
-1 . 3877e+003 
2 . 4228e+003 
1 . 9727e+003 
6 . 7734e+002 
8 . 6056e+002 
-1 . 8292e+002 


1. 5414e+001 
-1 . 2756e+002 
-2 . 7571e+002 
3 . 5902e+002 
9 . 1517e+002 
-1 . 1191e+003 
-1 . 3111e+003 
-2 . 5450e+002 
-2 . 5366e+003 
-4 . 9978e+003 
-6 . 3573e+002 
-2 . 0691e+003 
4 . 4304e+002 


-4 . 1640e+001 
-7 . 5931e+000 
1 . 4344e+002 
-2 . 6316e+002 
-1 . 9123e+002 
1. 6542e+003 
2 . 6818e+003 
-2 . 4940e+003 
3 . 7422e+003 
-1. 9697e+003 
2 . 8427e+003 
-2 . 5017e+003 
1 . 0769e+003 


-1 . 1395e+001 - 
-1 . 2589e+002 
-1 . 2242e+002 
7 . 4182e+001 - 
3 . 4736e+002 - 
3 . 6784e+002 
-6. 1243e+002 
-5 . 4361e+002 - 
-1 . 8956e+003 
-6 . 1339e+003 - 
-4 . 3431e+002 
-1 . 5847e+003 - 
1. 0643e+002 


Column 13 
-1 . 2417e+001 
2 , 1396e+001 
1 . 2171e+002 
-1 . 9983e+002 
-2 . 4419e+002 
2 . 0158e+002 
1 . 4490e+003 
-1 . 3702e+002 
4 . 8391e+003 
-7. 6598e+003 
5 . 0499e+003 
-1 . 4888e+004 
6 . 6320e+003 


1 . 3320e+002 
2 . 5421e+002 
3 . 3624e+002 
9 . 5962e+002 
1. 0894e+003 
2 . 2563e+003 
3 . 4127e+003 
2 . 6048e+003 
1 . 7062e+003 
9 . 3352e+002 
4 . 0328e+002 
5 . 9107e+002 
1 . 3 604e+002 


1 . 1287e+001 
2 . 5567e+001 
8 . 9461e+001 
1 . 609 le+002 
3 . 0721e+002 
3 . 9405e+002 
1 . 0956e+003 
9. 1853e+001 
3 . 4725e+003 
4 . 3988e+003 
2 . 4845e+003 
9 . 2850e+003 
4 . 8257e+003 
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Input Matrix [B] 

9 . 9293e+002 
-8 . 892 6e+002 
9 . 947 8e+002 
-7 . 0708e+002 
7 . 7543e+002 
-6 . 4590e+002 
3 . 3104e+002 
-2 . 7082e+002 
2 . 2690e+002 
1.4710e+002 
5 . 9450e+001 
8 . 4334e+001 
-1 . 8368e+001 


Output Matrix [C] 

Columns 1 through 6 

-6 . 8252e+002 -9.4337e+002 -4.1131e+002 8.2356e+001 6.7933e+002 

-7 . 2855e+002 1.4783e+002 8.9199e+002 -6.8364e+002 -4.6891e+002 


Columns 7 through 12 

9 . 3072e+001 -2.3150e+002 -1.4495e+002 -7.8913e+001 -1.3905e+002 
3 . 4574e+002 -1.6555e+002 2.3060e+002 -1.8832e+002 3.4148e+001 - 

Column 13 
3. 6269e+001 
-1 . 1410e+002 


Input/Output Matrix [D] 

0 

0 


4 . 0654e+002 
5 . 1561e+002 


2 . 2953e+001 
9 . 9862e+001 
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8.3.2 3D VDC Inlet Model: Downstream Mach Number Perturbation 


8.3.2. 1 15 x 3 x 3 Grid 


System Matrix [A] 


Columns 1 through 6 
8 . 9089e+001 


2 . 9149e+002 
■1 . 4482e+002 
3 .2395e+001 
■1 . 0587e+002 
3 . 9 694e+001 
4 . 3595e+001 
2 . 1687e+001 
1. 3773e+001 
1. 1264e+001 
5 . 3 664e+000 
6 . 7474e+000 
1. 6037e+000 
1 . 2031e+000 
4. 0208e-001 
3 . 3726e-001 
3 . 0765e-001 


-2 . 8477e+002 
-4 . 0097e+002 
8 . 2 662e+002 
-1 . 2348e+002 
4 . 1102e+002 
-1 . 5959e+002 
-1. 6895e+002 
8 . 4185e+001 
-5 . 3478e+001 
4 . 3730e+001 
-2 . 0851e+001 
2 . 6213e+001 
-6 . 2298e+000 
4 . 6735e+000 
1 . 5620e+000 
1 . 3101e+000 
-1 . 1951e+000 


-3 . 1694e+001 
-6 . 5305e+002 
-5 . 4169e+002 
1 . 9769e+002 
-9 . 4458e+002 
2. 8944e+002 
3 . 1327e+002 
-1 . 5718e+002 
1. 0076e+002 
-8 . 2295e+001 
3 . 9162e+001 
-4. 9248e+001 
1 . 1705e+001 
-8 . 7810e+000 
-2 . 9351e+000 
-2 . 4603e+000 
2 . 24 52e+000 


Columns 7 through 12 


1 . 9546e+001 
■2 . 7240e+001 
2 . 1635e+002 
1. 8545e+002 
■1 . 2541e+003 
1 . 0853e+003 
■2 . 1997e+003 
1 . 5388e+003 
1. 5335e+003 
1. 1549e+003 
5 . 4566e+002 
6 . 8607e+002 
1. 6192e+002 
1 . 2151e+002 
4 . 0631e+001 
3 . 4040e+001 
3 . 1072e+001 


8 . 4557e+000 
8 . 5437e+000 
-1. 1723e+002 
9 . 5964e+001 
5. 9512e+002 
-3 . 7373e+002 
1. 5421e+003 
-1 . 2908e+003 
2 . 1500e+003 
-1 . 4212e+003 
6 . 4794e+002 
-8 . 1534e+002 
1 . 9075e+002 
-1 . 4326e+002 
-4 . 7967e+001 
-3 . 9972e+001 
3 . 6597e+001 


-2 . 6019e+000 
-2 . 4834e+001 
-9 . 9812e+001 
5 . 0899e+001 
3 . 4117e+002 
3 . 3420e+001 
1 . 2776e+003 
-1 . 9337e+003 
-1. 9758e+003 
2 . 7745e+003 
-1. 6890e+003 
1. 9103e+003 
-4. 6009e+002 
3 . 4526e+002 
1 . 1557e+002 
9 . 6734e+001 
-8 . 8326e+001 


Columns 13 through 17 
1 . 5649e+000 


■5 . 5232e+000 
■1 . 7775e+000 
•4 . 1420e+000 
■1 . 4186e+001 
4 . 2598e+001 
•9 . 5518e+001 
1 . 0522e+002 
■3 . 0088e+000 
4 . 4558e+002 
9 . 7 643e+002 
3 . 2367e+003 
3 . 6401e+003 
4 . 0288e+003 
2 . 0885e+003 
1 . 5413e+003 
1 . 3061e+003 


1 . 4593e+000 
4. 8371e+000 
-3 . 6660e-001 
5 . 0086e+000 
2 . 0580e+001 
-3 . 9931e+001 
1. 1678e+002 
-1 . 4183e+002 
-7 . 9331e+001 
-3 . 4288e+002 
-8 . 3997e+002 
3 . 2056e+003 
5 . 9916e+003 
-8 . 3056e+003 
-7 . 5987e+003 
-5 . 7425e+003 
4 . 0209e+003 


-5 . 5831e-001 
-2 . 3537e+000 
-2 . 9280e+000 
-9 . 2137e-002 
3 . 9183e+000 
1 . 3587e+001 
-4 . 4895e+000 
6. 0799e+000 
-9. 5073e+001 
2 . 3974e+002 
3 . 57 64e+002 
-8 . 5378e+002 
-2 . 5077e+003 
5 . 3972e+003 
-1 . 8870e+003 
-7 . 0291e+003 
2 . 5996e+003 


-1 . 4626e+001 
-7 . 6740e-001 
9 . 9530e+001 
-7. 1802e+001 
1 . 0870e+003 
-2 . 1048e+002 
-1 . 8525e+002 
9 . 5383e+001 
-6 . 2517e+001 
5. 0883e+001 
-2 . 4234e+001 
3 . 0470e+001 
-7 . 2391e+000 
5 . 4311e+000 
1 . 8147e+000 
1 . 5216e+000 
-1 . 3885e+000 


9 . 0706e+000 
4 . 1180e+001 
6. 7762e+001 
-1 . 0573e+001 
-1. 3407e+002 
-2 . 2376e+002 
-2 . 6957e+002 
5 . 3048e+002 
8 . 0054e+002 
-2 . 3123e+003 
3 . 1305e+003 
-2 . 9995e+003 
6 . 5408e+002 
-4 . 9150e+002 
-1. 6386e+002 
-1 . 4433e+002 
1. 2768e+002 


-4 . 2748e-001 
-1. 5652e+000 
-7 . 8528e-001 
-9 . 2730e-001 
-2 . 5391e+000 
1 . 1078e+001 
-2 . 2706e+001 
2 . 9167e+001 
-1 . 2709e+001 
1 . 3096e+002 
2 . 5426e+002 
-8. 5345e+002 
-1. 9577e+003 
2 . 5718e+003 
2 . 6544e+003 
2 . 1499e+003 
-1 . 4804e+003 


-1 . 5213e+001 
2 . 1118e+001 
8 . 6629e+002 
-1. 0769e+003 
-1 . 4188e+003 
1 . 1077e+003 
1 . 2994e+003 
-6 . 0640e+002 
3 . 6381e+002 
-2 . 9995e+002 
1 . 4336e+002 
-1 . 8017e+002 
4 . 2844e+001 
-3 . 2139e+001 
-1 . 0739e+001 
-9 . 0137e+000 
8 . 2 194e+000 


5 . 3139e+000 
2 . 0660e+001 
1 . 8274e+001 
6 . 5951e + 000 
1 . 5380e+000 
-1.3792e+002 
1 . 7595e+002 
-1 . 8416e+002 
8 . 2780e+002 
-2 . 9372e+003 
-1 . 6625e+003 
3 . 3137e+003 
-1 . 0364e+003 
7 . 5270e+002 
2 . 43 3 8e+002 
2 . 1839e+002 
-1 . 9278e+002 


-1 . 4646e+000 
-4 . 0649e+000 
5 . 3497e+000 
-7 . 9436e+000 
-3 . 9295e+001 
4 . 1809e+001 
-1 . 8336e+002 
2 . 3211e + 002 
2 . 8491e+002 
1. 6351e+002 
7 . 5737e+002 
-3 . 9281e+003 
-7 . 4874e+003 
1 . 0405e+004 
1 . 1980e+004 
2 . 2744e+004 
-1 . 1146e+004 


3 

1 

3 

1 

-1 

-5 

-1 

6 

-3 

3 

-1 

1 

-4 

3 

1 

9 

-8 


-6 

-1 

5 

-2 

-1 

1 

-5 

6 
3 
1 
2 

-6 

3 

-2 . 
-7 . 
-6 
5 


. 9396e+001 
. 4638e+002 
. 6901e+001 
. 4102e+002 
. 0237e+002 
. 3801e+002 
. 4850e+003 
. 2848e+002 
. 7131e+002 
. 0688e+002 
. 4346e+002 
. 8097e+002 
. 3122e+001 
. 2343e+001 
. 0804e+001 
0699e+000 
2703e+000 


1248e+000 
9513e+001 
5233e+000 
3554e+001 
03 64e+002 
7364e+002 
4335e+002 
1092e+002 
4747e+002 
7109e+003 
139 6e + 003 
3908e+003 
3619e+003 
3743e+003 
6408e+002 
0228e+002 
7996e+002 
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Input Matrix [B] 

9.2761e+002 
-1 . 0753e+003 
7 . 7323e+002 
-1 . 6571e+002 
5 . 532 6e+002 
-2 . 0849e+002 
-2 . 2 685e+002 
1. 1292e+002 
-7 . 1737e+001 
5. 8665e+001 
-2 . 7953e+001 
3 . 5146e+001 
-8 . 3533e+000 
6 . 2 665e+000 
2 . 0943e+000 
1 . 7566e+000 
-1 . 6024e+000 


Output Matrix [C] 

Columns 1 through 6 

-6 . 2033e+002 -9.1893e+002 -7.1167e+002 6.5604e+001 3.4451e+002 

-6 . 8972e+002 -5.5828e+002 3.0226e+002 -1.5240e+002 -4.3299e+002 

Columns 7 through 12 

8 . 3295e+001 -4.8461e+001 -6.1499e+001 5.7470e+001 2.1684e+001 - 

-2 . 1101e+002 1 . 0260e+002 3.7078e+001 1.1810e+001 1.7691e+001 - 

Columns 13 through 17 

-4 . 3743e+000 2.9658e+000 -2.8631e+000 -1.3750e+000 -2.0617e-001 

-7 . 0344e+000 7.5784e+000 -1.3482e+000 -1.7679e+000 -1.0116e+001 


Input/Output Matrix [D] 

o 

0 


1 . 1812e+002 
1 . 7196e+002 


1 . 0069e+001 
3 . 3818e+001 
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8. 3.2.2 20 x 3 x 3 Grid 


System Matrix [A] 


Columns 1 through 6 
8 . 0831e+001 -2.3762e+002 


2 . 4335e+002 
1 . 2008e+002 
4 . 9785e+001 
8 . 9752e+001 
4 . 4009e+001 
3.5981e+001 
2 . 6242e+001 
7 . 8240e+000 
1 . 3911e+001 
4 . 1429e+000 
6 . 4412e+000 
1 . 5951e+000 
1. 1163e+000 
3 . 9209e-001 
5 . 0634e-001 
9 . 4549e-002 


-4 . 07 68e+002 
8 . 2492e+002 
-2 . 3172e+002 
3 . 9894e+002 
-2 . 0361e+002 
-1 . 6012e+002 
1 . 17 07e+002 
-3 . 4927e+001 
6 . 2035e+001 
-1 . 8490e+001 
2 . 8742e+001 
-7 . 1165e+000 
4 . 9805e+000 
-1 . 7494e+000 
2 . 2591e+000 
-4 . 2184e-001 


- 2 . 4279e+001 
-6 . 5485e+002 
-4 . 3398e+002 
2 . 6493e+002 
-7 . 7778e+002 
3 . 2216e+002 
2 . 5023e+002 
-1 . 8412e+002 
5 . 5503e+001 
-9 . 8554e+001 
2 . 93 08e+001 
-4 . 5573e+001 
1 . 12 84e+001 
-7 . 8977e+000 
2 . 7738e+000 
-3 . 5820e+000 
6 . 6911e-001 


-1 . 1006e+001 
6 . 5546e+001 
1 . 9499e+002 
-1 . 6325e+002 
1 . 0881e+003 
-3 . 4777e+002 
-2 . 2643e+002 
1 . 6992e+002 
-5 . 2264e+001 
9 . 2402e+001 
-2 . 7469e+001 
4 . 2704e+001 
-1 . 0573e+001 
7 . 3986e+000 
-2 . 5990e+000 
3 . 3557e+000 
-6 . 2 672e-001 


-9 . 3918e+000 
3 . 0660e+001 
7 . 2874e+002 - 
-1 . 0878e+003 
-1 . 2054e+003 
1 . 0986e+003 - 
1 . 0898e+003 - 
-7 . 4159e+002 
2 . 063 6e+002 - 
-3 . 7090e+002 
1 . 1095e+002 - 
-1 . 7250e+002 
4 . 2728e+001 - 
-2 . 9905e+001 
1 . 0505e+001 - 
-1 . 3564e+001 
2.5326G+000 - 


Columns 7 through 12 


1. 6391e+001 
2 . 8012e+001 
1 . 7546e+002 
2 . 1594e+002 

1 . 03 92e+003 
1 . 2189e+003 
1 . 5809e+003 
1. 6258e+003 
8. 0071e+002 
1 . 2549e+003 
3 . 6853e+002 
5 . 6952e+002 
1 . 4027e+002 
9 . 8165e+001 
3 . 4460e+001 
4 . 4519e+001 
8 . 3254e+000 


7 . 8957e+000 
1. 6930e+000 
-1. 5057e+002 
1 . 6919e+002 
7 . 3341e+002 
-6 . 4558e+002 
1 . 5800e+003 
-1 . 9511e+003 
1 . 7172e+003 
-2 . 2654e+003 
6 . 2851e+002 
-9 . 6798e+002 
2 . 3737e+002 
-1 . 6595e+002 
5 . 82 57e+001 
-7 . 52 60e+001 
1 . 4078e+001 


3 . 5453e+000 

6 . 3 881e+000 
-3 . 9959e+001 

5 . 03 04e+001 
1 . 9522e+002 

-1 . 7947e+002 
8 . 0407e+002 
-1 . 7080e+003 
-5 . 3548e+002 
2 . 0025e+003 
-5. 3716e+002 
8 . 4994e+002 
-2 . 1770e+002 
1 . 5333e+002 
-5 . 3606e+001 
6. 9220e+001 
-1 . 3171e+001 


9 . 4518e+000 
5 . 3172e+001 
8 . 9316e+001 
-5. 6047e+001 
-2 . 3762e+002 
-1 . 013 7e+002 
-5 . 053 0e+002 
1 . 3627e+003 
-2 . 2080e+002 
-3.797 le+003 
3 . 1795e+003 
-3 . 7259e+003 
8 . 5394e+002 
-6 . 0817e+002 
2 . 1403e+002 
-2 . 7722e+002 
5 . 0869e+001 


4 . 1806e+000 
1 . 7608e+001 
7 . 6169e+000 
6 . 0994e+000 
1 . 8842e+001 
-1 . 0716e+002 
1 . 8952e+002 
-2 . 4338e+002 
-1. 9028e+002 
-2 . 5941e+003 
-9 . 5163e+002 
2 . 2666e+003 
-7 . 7236e+002 
5. 1506e+002 
-1 . 8322e+002 
2 . 3313e+002 
-4 . 2669e+001 


Columns 13 through 17 
1 . 2932e+000 


6 . 6970e+000 
9 . 0156e+000 
4 . 5313e+000 
2 . 0322e+001 
1 . 9609e+001 
2 . 1608e+001 
7 . 5223e+001 
3 . 1195e+001 
8 . 3095e+002 
6 . 2929e+002 
2 . 6414e+003 
2 . 3172e+003 
3 . 2785e+003 

1 . 03 64e+003 
1 . 4858e+003 
3 . 0408e+002 


8 . 8110e-001 
2 . 9214e+000 
-2 . 7423e+000 
5 . 4978e+000 
2 . 1075e+001 
-3 . 1577e+001 
8 . 8711e+001 
-1 . 3838e+002 
-1 . 4075e+002 
-4 . 7256e+001 
-4 . 3347e+002 
2 . 4026e+003 
4 . 1506e+002 
-7 . 3598e+003 
-5 . 7799e+002 
-6 . 2448e+003 
-5 . 2502e+001 


-1. 4053 e+ 000 
-5 . 0179e+000 
2 . 6890e+000 
-6 . 9861e+000 
-2 . 6298e+001 
4 . 5523e+001 
-1 . 2143e+002 
1 . 9599e+002 
2 . 0121e+002 
1 . 7463e+002 
6 . 9859e+002 
-3 . 7595e+003 
-1 . 7694e+003 
1 . 8380e+004 
-2 . 5993e+003 
7 . 5030e+003 
~1.5795e+003 


2 . 3677e-001 
1 . 7242e-001 
-4 . 1457e+000 
4 . 729 le+000 
1. 8847e+001 
-1. 5061e+001 
6 . 3272e+001 
-1 . 133 5e+002 
-1 . 0326e+002 
1 . 9110e+002 
-1. 3583e+002 
1 . 0895e+003 
-6 . 6328e+002 
-6 . 2113e+003 
-4 . 3101e+002 
-1 . 1786e+004 
-2 . 653 6e+003 


-6 . 2074e-001 
-2 . 0624e+000 
2 . 0549e+000 
-3 . 9265e+000 
-1 . 4973e+001 
2 . 1823e+001 
-6 . 3606e+001 
1 . 0600e+002 
1 . 0581e+002 
2 . 5657e+001 
3 . 1492e+002 
-1 . 7664e+003 
-5 . 82 16e+002 
8 . 9025e+003 
-3 . 3418e+002 
8 . 6370e+003 
2 . 8574e+002 


3 . 8519e+001 
1 . 4726e+002 
9 . 0741e+001 
2 . 5794e+002 
1. 9803e+002 
6 . 3949e+002 
1 . 3591e+003 
8 . 2990e+002 
2 . 1880e+002 
4 . 0243e+002 
1 . 1872e+002 
1 . 8520e+002 
4 . 5956e+001 
3 . 2154e+001 
1. 1297e+001 
1 . 4585e+001 
2 . 7223e+000 


6 . 0721e+000 
2 . 3 139e+001 
2 . 1882e+000 
2 . 1821e+001 
8. 0148e+001 
1 . 8433e+002 
4 . 2139e+002 
6 . 0204e+002 
6 . 0454e+002 
1 . 7795e+003 
2 . 0955e+003 
6 . 9498e+003 
3 . 9212e+003 
2 . 4785e+003 
8. 8212e+002 
1 . 0918e+003 
2 . 12 44e+002 
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Input Matrix [B] 

8 . 6027e+002 
-9 . 7457e+002 
6 . 5451e+002 
-2 . 6276e+002 
4 . 7994e+002 
-2 . 3 679e+002 
-1. 9174e+002 
1 . 3995e+002 
-4 . 1738e+001 
7 . 4189e+001 
-2 . 2095e+001 
3 . 4352e+001 
-8 . 5060e+000 
5 . 9528e+000 
-2 . 0909e+000 
2 . 7 OOle+OOO 
-5 . 0419e-001 


Output Matrix [C] 

Columns 1 through 6 

-5 . 9238e+002 -8.4183e+002 -5.8859e+002 1.5571e+002 3.0401e+002 

-6 . 2388e+002 -4.9035e+002 2.8578e+002 -2.1171e+002 -3.7135e+002 

Columns 7 through 12 

6 . 4280e+001 -6.8326e+001 -1.4187e+001 7.4199e+001 1.4560e+001 - 

-1 . 8064e+002 1.2258e+002 3.9477e+001 -1.0981e+000 1.6838e+001 - 

Columns 13 through 17 

-8 . 3285e+000 5.8599e-001 -1.9637e+000 -1.7740e+000 -3.7570e-001 

-1 . 5866e+000 5.9150e+000 -8 . 4767e+000 3.4320e+000 -4.2146e+000 


Input/Output Matrix [D] 

o 

0 


5 . 9407e+001 
2 . 2936e+002 


1. 3534e+001 
3 . 1703e+001 
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8.3.3 ID Model from 2D Averaged Data 


System Matrix [A] 


Columns 1 through 6 
■9 . 8584e+001 


2 . 8000e+002 
•3 . 03 12e+002 
■2 . 4737e+002 
■8 . 8897e+001 
1 . 5402e+002 
■1.6297e+002 
1 . 6455e+002 
1.0043e+002 
1.6757e+001 
8 . 3077e+001 
1 . 0901e+000 
4 . 8729e+001 


2 . 4455e+002 
-1 . 4186e+002 
-1 . 5407e+003 
-2 . 4315e+002 
-2 . 2576e+002 
3 . 8103e+002 
-2 . 3407e+002 
2 . 9834e+002 
-1. 6631e+002 
2 . 6745e+001 
-1 . 4773e+002 
-5 . 9525e-001 
-8 . 4991e+001 


Columns 7 through 12 


1 . 2185e+002 
-2 . 3521e+002 
4. 6027e+002 
-5 . 6916e+002 
2 . 8608e+002 
-3 . 2350e+003 
-8 . 7731e+002 
3 . 4071e+003 
-8 . 9785e+002 
1 . 3322e+002 
-1 . 0200e+003 
-1 . 255'6e+001 
-5. 6226e+002 


-4 . 4489e+001 
-2 . 0590e+002 
5 . 0430e+002 
-6 . 0448e+002 
-4 . 4595e+002 
8 . 7985e+002 
-3 . 0424e+003 
-1 . 9011e+003 
4 . 5798e+003 
-6 . 0813e+002 
2 . 1877e+003 
-4 . 6541e+001 
1 . 0281e+003 


Column 13 
4 . 9590e+001 
-6 . 2108e+001 
-1. 1102e+002 
-1 . 9024e+002 
1 . 3820e+002 
-2 . 7672e+002 
-5 . 0154e+002 
7 . 0727e+002 
3 . 5872e+002 
8 . 7032e+002 
4 . 3775e+003 
-6 . 5665e+003 
-2 . 0530e+003 


1 . 6793e+002 
1 . 4511e+003 
-3 . 9981e+002 
-1 . 3903e+003 
-6 . 2326e+002 
6 . 8001e+002 
-5 . 9834e+002 
4 . 9248e+002 
-3 . 4549e+002 
5 . 4938e+001 
-2 . 77 94e+002 
-2 . 6441e+000 
-1. 6060e+002 


2 . 7052e+001 
-1. 5972e+002 
3 . 1377e+002 
-4 . 0688e+002 
-3 . 5490e+001 
-3 . 6715e+002 
-9 . 0788e+002 
-4 . 4534e+003 
-1 . 4257e+003 
1.703 le+002 
-3 . 3635e+003 
-1 . 4002e+002 
-1 . 7322e+003 


2 . 0341e+002 
-2 . 4207e+002 
1 . 2190e+003 
-4 . 8376e+002 
-2 . 2527e+003 
1 . 6473e+003 
-5 . 5400e+002 
8 . 7442e+002 
-4 . 4008e+002 
6 . 8069e+001 
-4 . 1020e+002 
-3 . 7 015e+000 
-2 . 3313e+002 


-2 . 8182e+001 
4 . 8691e+001 
-4 . 2673e+001 
1 . 4473e+002 
-5 . 13 68e+001 
3 . 2932e+002 
3 . 0846e+002 
6 . 5000e+002 
4 . 42 04e+002 
-3 . 4059e+001 
3 . 9685e+003 
• 1 . 7375e+002 
7 . 3388e+002 


-7 . 1015e+001 
1 . 2358e+002 
5 . 6456e+002 
2 . 0570e+003 
-1 . 4243e+002 
3 . 4002e+002 
-5 . 3727e+002 
5 . 3415e+002 
-2 . 7518e+002 
4. 7297e+001 
-2 . 3463e+002 
-1. 6848e-001 
-1 . 2837e+002 


-8. 6517e+001 
1. 0008e+002 
2 . 7896e+001 
3 . 0481e+002 
-1. 9536e+002 
8 . 4099e+002 
6 . 7582e+002 
-6 . 2777e+002 
2 . 5041e+003 
-4 . 4839e+003 
-2 . 1793e+003 
-2 . 1985e+002 
-3 . 9973e+003 


1 . 7131e+002 
3 . 4936e+002 
3 . 5812e+002 
1 . 6233e+003 

2 . 3 804e+002 
6 . 8877e+002 

3 . 2 588e+003 
1 . 5768e+003 
8 . 5709e+002 
1 . 6234e+002 
6 . 1397e+002 

1 . 3 009e+001 
3 . 8545e+002 


2 . 1196e+001 
1 . 7938e+001 
1 . 2451e+002 
7 . 1453e+001 
8 . 0668e+001 
2 . 543 6e+002 
1. 5670e+002 
9 . 0295e+002 
5 . 0311e+002 

2 . 3 094e+002 
1 . 0737e+003 
1 . 4654e+002 
3 . 1347e+003 
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Input Matrix [B] 

9 . 7210e+002 
8 . 053 6e+002 
1. 1889 e+ 003 
1.0025e+003 

5 . 03 62e+002 
-8 . 5236e+002 

7 . 6198e+002 
-8 . 0286e+002 
4 . 9 182e+002 
-8 . 0081e+001 
4 . 1501e+002 
3 . 5179e+000 

2 . 3 900e+002 


c 


Output Matrix [C] 


Columns 1 through 6 
-6.3659e+002 8.1445e+002 -7.0900e+002 

-7.3606e+002 6.3800e+000 9.6624e+002 


1 . 0345e+003 -7.9062e+001 
6 . 2264e+001 -5.0451e+002 


Columns 7 through 12 
8 . 0468e+002 4.8980e+002 

-3 . 4003e-002 -6.5999e+002 


4.5218e+002 -1.6378e+002 -3.0209e+002 - 
-2 ! 3037e+002 -3.8097e+001 -2.9737e+002 


Column 13 
1 . 2770e+002 
2 . 3890e+002 


'n 


Input/Output Matrix [D] 

o 

0 


6 . 4761e+002 
6 . 9652e+002 


6 . 4401e+001 
1 . 9047e+002 
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